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SUMMARY 


The general airfoil theory is extended to cover the case of air- 
foils having air inlets and suction slots. The analysis, employing 
formal potential theory, rests on the usual assumptions of an in- 
compressible, nonviscous fluid satisfying the Kutta condition. 
Subject to these assumptions, it is found that the lift and drag of 


such an airfoil are given by 


Bre 
Cr = 2r(a te) +2 Co» tan (: es ) 
n - 
Cp = 2 pm Con 
n 


where the symbols are defined on Fig. 6. 

These results have been confirmed by simple tests, the re- 
sults of which are summarized in Fig. 8. 

In general, it is indicated that an air inlet or suction slot on 
the upper surface, +, increases the lift or, on the lower surface, 
—8, decreases the lift in a manner analogous to changing the 
camber of the airfoil. An air inlet near the nose, 8 = 0, is indi- 
cated to have no effect on the lift, while a suction slot on the upper 
surface, or the equivalent of a source on the lower surface near 
the trailing edge, is indicated to give large lift increments. The 
latter arrangements are considered to be of possible interest as a 
control or high lift device. All of these effects are hydrodynamic 
and independent of the effect of the inlet or slot on the boundary 
layer. 


NOTATION 


z= Cartesian coordinates 

free-stream velocity 
component of velocity along x-axis 

= component of velocity along y-axis 
component of velocity along s-axis 
general displacement s = V x? + 9' + s* 
general velocity g = Vu? + 2° + w* 
polar coordinates (see Fig. 3) 
radius of cylinder 

= chord of airfoil 
polar angle of source-sink discontinuity (see Fig. 1) 
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e + 7 = polar angle of point transforming into trail- 
ing edge of airfoil 

angle of attack 
velocity volential 
stream function 

= circulation 
complex coordinate (see Fig. 1) 
complex coordinate of poles of Eq. (10 
potential function 
complex coordinate (see Fig. 1) 
complex operator 
coordinate of source in z-plane (see Fig. 1) 

= volume flow from source (cu.ft. per ft. per sec.) 

= m/ Uc, flow coefficient (positive for sink) 

= absolute pressure 

= density 

2 = contours of integration (see Fig. 2) 
», 53, Se = surface of integration (see Fig. 2) 
= net aerodynamic force on airfoil in chordwise direc- 
tion 

net aerodynamic force on airfoil normal to chord line 
constants defined by Eqs. (10) and (11) 


(1) INTRODUCTION 


ie 1912 JOUKOWSKY PUBLISHED a treatise! presenting 
for the first time a complete and practical solution 











Flow in {-plane obtained from flow in z-plane by con- 


Fic. 1. 
formal transformation defined by Eq. (1). 
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for the lift and drag of an airfoil. This work, conclud- 
ing a long quest for such a solution, was partly based 
on, and partly done concurrently with, that of Kutta 
and Lanchester. Although refinements and extensions 
have since been made (notably by Glauert,? Theodor- 
sen,* and Allen‘), the methods and assumptions pre- 
sented in this original paper constitute the basis of air- 
foil theory as applied today. 

This theory has long served as a rigorous framework 
around which specific experimental investigations 
could be organized and which, in recent years, has been 
successfully used as an analytical tool in the design of 
high-performance sections. 

In the present paper the general theory is extended to 
cover the case of the airfoil having air inlets and suction 
slots on its boundary. This solution is obtained by 
representing the suction slot or air inlet by one or more 
sink discontinuities. The formal problem herein con- 
sidered thus consists of solving for the lift and drag of an 
airfoil having one or more source-sink discontinuities 
on its boundary. 

The reader is referred to references 2 and 5 for the 
necessary background for the development that fol- 
lows. 


(II) FORMAL SOLUTION 


(A) Potential Function for Cylinder with Source-Sink 
Discontinuity on Boundary 


The transformation 
g = —iR[(z + 1)/(2 — 1)) (1) 


takes the imaginary axis, zy, in the z-plane into a circle 
of radius Rin the {-plane. (See Fig. 1.) In particular, 
the points —i, 0 and +7 are transformed into points at 
R, + iR, and —R. If, therefore, a simple flow can be 
built up in the z-plane in which 7y is a streamline and in 
which there is a source-sink discontinuity located on zy, 
then transforming this flow with Eq. (1) will yield the 
flow about a cylinder with a source or sink on its bound- 
ary. 

If a source-sink discontinuity is simply placed on the 
iy-axis, say at + /1, zy will be a stream line which will, as 
desired, become a circular cylinder. However, if a 
single source or sink is located in the z-plane, a singular- 
ity of opposite sign is automatically placed at infinity 
(continuity requirement). Since the point at infinity 
in the z-plane transforms into a point on the cylinder 
in the ¢-plane, the transform of a single sink at +/7 is a 
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cylinder with a sink and a source on its boundary. To 
overcome this difficulty, assuming a sink at /7, a source is 
placed at +1, and, in order to preserve zy asastream line, 
a second image source is placed at —1. The strength of 
the sources at —1 and +1 are noted as m and that of 
the sink is, on this basis, —2m.* With this arrange- 
ment there is no singularity at infinity. Thus, in the 
¢-plane, a cylinder is produced with a sink of strength 
— 2m on its boundary, a source at the origin (the trans- 
form of z = —1), and a source at infinity (transform of 
z= +1). 

The potential function of the new flow could be 
found by substituting the inverse of the transformation 
function (1) into the potential function of the simple 
source-sink flow in the z-plane. In this case, however, 
the new flow is composed of point sources and sinks of 
known strength and location, and it is simpler to write 
down the potential function directly. Thus since the 
potential function of a source (+) or a sink (—) located 
at the origin and of strength m is 


w= +(m/2zx) Ing (2) 
that of a sink of strength —2m located at Re” is 
w = —2(m/2m) In (2 — Re'®) 


The total potential function for the flow in the £-plane 
is thus 


m , m 
oe “82h g = Re*®) + 5, ns 
ar aT 


A uniform flow of velocity Up from the right may be 
added by superimposing the potential function for 
uniform flow about a cylinder.! Since the cylinder 
stream line is common to both flows, this boundary is 
unaffected by the addition. The total potential func- 
tion for uniform flow past a cylinder with a sink of 
strength —2m in its boundary is 


R 
a ~ [2 In (¢§ — Re®) — nt] — Up (: + ‘Yo 
2r f 


This can be transformed into the flow about such a 
cylinder at angle of attack a by the transformation 
3' = te—®, or the new potential function can be ob- 
tained by substituting the inverse, ¢ = 2’ '* into Eq. 
(3). That is, 





* Note, however, that the volume removed from the flow ex- 
ternal to the cylinder is simply m. 





m ; x 
w= — —[2 In (s’* — Ref) — In 2’e*]— Us (s"* + ra) “ 
anv ‘ 


” 


where the a is added to the Re’ + ™ term in order to keep the sink at the same point on the cylinder, irrespective 


of the angle of attack a. 


Finally, the effect of a circulation about the cylinder can be included by superimposing the potential function 


for a point vortex at the origin. Eq. (4) thus becomes 
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m , ; ‘ R? 4 ; 
w= — [2 In (2'* — Re“®+@)) — In ge] — Up ( se + —) — i — In ze (5) 
2a ze* 2a 


where the primes have now been dropped from the s’s for simplicity. 


Adding more source-sink discontinuities merely adds additional terms of the form of the first of Eq. (5). 


fore, for the general case, 


w= 


(B) Lift and Drag of the Cylinder with One or More 
Source Sink Discontinuities on Its Boundary 
The pressure forces acting on any body whose sur- 
face is a stream line and for which the potential func- 
tion for the flow is known can be computed directly by 
application of the Blasius formula, 


X-iY=(1, 2)pi J (dw/ds)*ds (7) 


In the present case, there being one or more source- 
sink discontinuities on the surface of the body, this 
boundary is not a stream line and Eq. (7) is not directly 
applicable. If, however, it be recalled that the sources 
and sinks are themselves being used to represent slots 
of finite width and if the momentum balance be written 
for this actual system, this difficulty can be overcome 
and a practical significance given to the solution of 
Eq. (7). Let the system be defined as that bounded 
by surfaces S;, Se, S3, and Sy (see Fig. 2), where S; is the 
exposed surface of the body including the walls of the 
slot or slots, Sz; is the exposed bounding surface, and 
5S; and S; are the side forces of the bounding surface. 
The x- and y-components of force are then 


~ on ft dy; —iY + if 7" 


where X and ¥ are the pressure forces on the body sur- 
face S; including those acting on the inside of the slot. 
Since S; is now a stream line, the ratio of momentum 


transfers are 
f (ypu dy — upv dx) 
€ 
2 


if (vpu dy — vpv dx) 
9 


Equating the forces to the rates of momentum change 
and simplifying, one obtains 


+ ) ‘ - a [2 In (se'* — Rest) — In ze] — Uy (c's + 
Lae av 


There- 


R? 


4 
2e 


) — In zal (6) 


2a 














System used in applying Blasius equation to body hav- 
ing source-sink discontinuity on boundary. 


X-1tY= (1/2)p1 D (dw/dz)*dz 


where X and Y are the pressure forces not only on the 
exposed surface of the body but also the walls of the 
slot and w is the potential function for the slot flow. 
The width of the slot (or slots) may, however, be re- 
duced to the point where the potential function w for 
the slot flow is, along the boundary C:, indistinguishable 
from the potential function (5) for the line sink. Thus 
Eq. (5) may be substituted into Eq. (7). In general, 
if a slot is represented by a line sink, the solution of the 
integral (1/2) 6 (dw/dz)* dz yields a quantity the 
“2 

real and imaginary parts of which are, respectively, the 
X- and Y-components of force on the body,- including 
the internal forces required to divert the internal flow 
parallel to the spanwise axis. 


Fic. 2. 





dw m ( 2 ") 
dz Qn \s — Re® 


() Ss 
dz} \ 


Expanding and rearranging, 


. R? :, : 5 
— Ure (1 -— “" —4 (8) 
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(=) - Uve™™ + | - 2 (=) Ure + 2 Ure iF + | (=) - 2Ure™*R? — 2 () + | + 
dz 2n 2r z 2r 2ar\2r 4e? Jj3° 
2 () UeltR? — 2% Vier: | 4 (Uy?e~ **R*) = 4 (5) vi _— x + 
Qa 23 2° s* 2r (sc — Re 
m iT (/m l m on | m \* l 
BS +s ae Re®) * F< + (Fe) Ue | sere +45) aoe 
The integration of Eq. (10) according to Eq. (7) involves the termwise integration of the following functions of 
K,; Ke Kz Kz Ks Ke Ky Ks 
9 meen mma ’ ——» —__————_) ——— -9 ; (11) 
z 2? 23 gt (g — 2) (f — 2)? 2(8 — 20) 2°(3 — 2p) 
one obtains 
K K l Ys 2 2"! : 
o"(5 — 2 ‘leow-i- 27? 7. =| 
(14 
from which it is apparent that (13) is itself composed of 
a series of terms of the form of (12). Thus, all of the 
functions of s appearing in Eq. (10) are, or can be, 
written in the form A/(z — 2)" 
Eq. (7) requires that these terms be integrated 
around a path C; enclosing poles (divisions by zero) at z 
= 2 = Oandz = % = Re*. Now the integral of each 
function around C, can be written as the integral around | 
Fic. 3. , f(z)dz C; plus the integrals around C, and C;, where C, and 
C; enclose the poles (see Fig. 3). The sum of these 
It is immediately apparent that all of these terms are of — line integrals, C; + Cy + C;, equals that around G 
one of two general types But path C; is indented so as to enclose no poles and, 
K/(¢ — )" (12 hence, by Cauchy's integral theorem is zeTO. The in- 
K/2"(s — x) (13) tegrals of tS) along paths C, and Cs enclosing the poles 
at 0 and Re® are obtained in each case as follows. Let 
Furthermore, breaking (13) up into rational fractions, 2 — 2 = re” (see Fig. 3.). Then for  # | 


dz ; “ve id ; De atte of ; 
a= K -~ = K 0 (15) 
° 0 yn tn a ] 
= re L? l(l—n 


and for = 1 


' dz : of ,' 
Kf. ~ = f, 1d6 = A2zm 16) 
\3 = Zo) 


Thus, by Eq. (15) the integrals of all terms for which m ~ 1 around C, and C; are zero, and by Eq. (16) the integral 
of all terms for which m = 1 is equal to K27i about the path enclosing its pole—i.e., about C,or C;. Since the in- 
tegral around C; equals that around C; + C, + C; and since that around C; is always zero, it follows that the i1- 
tegral of all terms around C; for which n + 1 is zero—i.e., 


G.IK dz/(z — 2)" = 0, ¥ 1 li 


and that of those for which m = 1 is A2ri—.e., 


GIK dz/(z — 2)] = K2zi (18 


Substituting Eq. (18) into the integral of Eq. (14), it follows that 
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K dz K ‘ A 
G., a = ri — 2ri - 0-0... — 0) = 0, forn 21 (19) 
C1 2"(g — 2) 0 
Using Eqs. (17), (18), and (19) it follows that the integral around C; of all terms of Eq. (10) are zero except the 
second and sixth. Thus, the integration of Eq. (10) reduces to 


dw\? a oe, 2 dz We. 2 dz 
Sg. (=) dz = | -2(%)e e+ 3. Ue |%. . + [4(™)c oe | Cs (g — Re’) 





or 


aii - 


dw\? w\.. i _ Oe A m\. ; Per a = 
gf. ( ~) dz =| -—2 Ue’* + 2—Upe™ |2ri + | 41 — Je oe” |2at = ZWnU oes — 2T Uce™ 
C: \ dz 2r 2a 27 


and Eq. (7) reads 


= ; dw\* ; _ : lee - ; 
X -1Y = (1/2) f (*) dz = (1/2)p1(2m Uoe’*t — 2TUce’*) = —pmUve'™™ — pr Ucei = — pUom X 
az 
(cos a + 7sin a) — ipl’ Up (cos a + 7 sin a) = (—pUom cos a + pl Up sin a) — 


i (pUo m sin a + pI' Up cos a) 
or 


X = —plUym cos a + pl'Uo sin a | 
, y S | 77 { (20) 
Y = pUym sin a + pI'Up cos a 
Eqs. (20) indicate that for the cylinder the components of force normal and parallel to the free stream direction—.e., 
the lift and drag are, respectively, pU) T and pUgm. (See Fig. 4.) 
Examining Eq. (S) in retrospect, it is apparent that in squaring dw/dz the only terms that do not produce zero 


forms are the products with — L'ye"*, and in squaring each of these occurs twice. Thus, Eq. (9) could be written 
(dw\° aie m 2 ] ae : ; y . 
( —} = —2U re - —— — t— + (terms which on integration = zero) (9a) 

° 9 s—R ip - 9 “ 

dz au Z é s om Z 


If now two or more source-sink discontinuities are located on the boundary, the potential function corresponding 
to Eq. (5) is as noted by Eq. (6). The derivative corresponding to Eq. (8) is 


dw My 2 ] aa te = . 
~=>|- m( —* = - >)| - ve (1 - 3 ° )-i5: (8a) 


As in the case of Eq. (9a), the square of Eq. (Sa) produces all zero forms except for the products with — Use"*,* 





which again each occur twice. Thus, 


leo 
atin 


dw\? : m 2 1 Sa ; : ; 
( ) = —2l | E( —~ +--- moi 6 + (terms which on integration = zero) (9b) 
dz * 2rz — Re z 


* The only form not previously discussed which develops is in the cross products between the source-sink potential functions. 





Here one obtains products of the form ¢ [dz/(z — a)(z — 8)] (a + B inside c). By breaking this up into rational functions, 
however, one obtains the difference of two integrals of the form of (12). This form is therefore also zero. 
which when substituted into Eq. (7) yields iy 


X -—iyY = (iI 2)pi(2Uvie'*>-m, — 21 Ue’) 


n 


or 





Lift = pUST | 
Drag = pUy) om, 
n 








(21) 





| 
| 
| 





From these results the following conclusions can be 
drawn: (1) the net forces on a body having one or more 
Source-sink discontinuities on its boundary are inde- Fic. 4. Showing resultant forces defined by Eq. (20). 
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Fic. 5. Showing how circulation is increased by suction slot. 


terminate; but, if it be assumed that the fluid taken 
in or discharged from the body flows initially or finally 
along a spanwise path, then the external plus the inter- 
nal forces are finite and given by Eq. (21). In particu- 
lar, (2) the lift of a cylinder with one or more source- 
sink discontinuities on its boundary is, as in the case of a 
cylinder without such discontinuities, given by L = 
pUl—i.e., the discontinuities themselves have no di- 
rect effect on lift, the external pressure changes being 
canceled by the internal pressure or inertiaforces. (In 
Section (D) it will be shown, however, that for airfoils 
the Kutta condition couples the source strength to the 
circulation and thus indirectly results in a change in 
lift). (3) The drag of a cylinder with one or more 
source-sink discontinuities on its boundary is equal to 
the rate of change of momentum in the free stream. 
(It should be noted that this drag will be zero if the 
air taken in or discharged is subsequently discharged or 
previously taken in from the free stream.) 

(C) Applicability of Results to an Airfoil 


A circle may be transformed into any airfoil by 
means of a transformation function of the form 


a2 a3 Qn (22 


a1 
jvet ses. ...2 


The best known systematic method for carrying out 
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such a transformation is due to Theodorsen and is de. 
scribed in reference 3. Having the potential function 
w of the flow about a cylinder with a source-sink dis. 
continuity on its boundary—i.e., Eq. (5)—one may 
determine the potential function for the flow about any 
airfoil, obtained by use of transformation (22), by sub. 
stituting into Eq. (5) the inverse of the particular 
transformation used. The lift and drag of the airfoil 
may then be found using Eq. (7), 


2)pi Ph (dw/ds)*at 


It is simpler, however, to obtain the required (dw/d¢) 
directly by multiplying the first derivative of the po- 
tential function for the flow about the cylinder, (dw/dz) 
by the reciprocal of the derivative of the transforma- 
tion function (dz/dt)—i.e., 


xX — 1Y = (23) 


(dw/d¢) = (dw/dz)(dz/dé) (24 


The lift and drag of the airfoil are then given by 


X -iY = = 1/21 B, pro dg = 























Now from Eq. (22) 





de 5 as 2a Bay 
dz 2? 3 z4 stdin 
or 
dz/d¢ = 1 + (a;/2? 


Substituting into Eq. (25) 


‘ ae dw\?/ a 

X -iY =(1 2) pi D ( “) (1 +—.. ) dz (26 
ce \dz 2? 

Substituting for (dw/dz)? from Eq. (9b), one finds that 

all integrals involving products with the a/z” terms 


are zero, so that Eq. (26) reduces to 


= (1/2)pi if. (dw/dz)*dz 


But the integral on the right is that carried out for the 
cylinder. Thus, the forces acting on an airfoil are th 
same as those acting on a cylinder having the same circ 
lation T and source strength m. The conclusions, stated 
previously, of the effect of a suction slot on the forces 
on a cylinder are fully applicable to the airfoil—namely, 
a source-sink discontinuity of itself produces no lift but 
if, through some means it can be made to change the circ 
lation about the airfoil, a change in lift will indirectly bt 
produced. The effect of the suction slot on the circul 
tion is considered in the next section. 
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(D) Effect of a Source-Sink Discontinuity on Circulation 


Flow around a sharp corner theoretically results in 
an infinite velocity. In the presence of a boundary 
layer through which pressure waves can propagate, 
an infinite velocity (and the associated pressure gradi- 
ents) cannot be produced as a steady state. The only 
alternative condition that can be established is for a 
stagnation point to occur at the trailingedge. That the 
latter condition actually exists was first postulated by 
Kutta and has subsequently been experimentally 
demonstrated to correspond to reality. An airfoil 
develops lift as a result of the development of sufficient 
circulation to so place the stagnation point on the trail- 
ing edge. Therefore, to determine the lift on an air 
foil, one calculates the circulation required to place the 
stagnation point on the trailing edge and substitutes 
this value into the general solution L = pU,I' [see 
Eq. (21)]. 

Starting with Eq. (8), one may proceed to determine 
the circulation required to place the stagnation point 
at the trailing edge as follows: 


dw m ( 2 *) U; (1 i, “] 
- a a < ————— — - _ é — -@e ~— 
dz In\z — Re® . 2? 


This must be zero at the point on the cylinder which 
transforms into the trailing edge; or, it must be zero 
for 
g= Ri@ +o Pe — Re** 
where (7 + €) is the value of 6 for the particular point 
that transforms into the trailing edge. (e€ is zero for 
symmetrical airfoils.) 
Substituting into Eq. (8b) 


dw m ( 2 4 l Unei® se 
-=— a 9 — — e 
dz 2a\— Re — Re® ~ Re** : 
Re _.. ir 1 ao 
( ~ Pogrit © ) Ree 9 (27) 


and solving for I’, one obtains 


f 9) 1é 

; e <e ' = 

rT = —im( - rc iB ote 1) ~i2eRL » puttin x 
=F 


( 1 —e ila +) 


first term by 


(28) 
Dividing top and bottom of the 


e‘[(e + 8)/2] and the second by 2%, 


r= —m tan(“ = *) + 4rRU, sin (a + €) (29) 


~ 


Similarly, for an arbitrary number of discontinuities 





— B 
r= }— m, tan (‘ — ) + 4rRUp sin (a + €) | (29a) 
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The second term in Eqs. (29) and (29a) may be 
recognized as the value of I’ for a normal airfoil having 


no discontinuities ori its surface. It thus appears that, 
af the circulation is determined by the Kutta condition, a 
source-sink discontinuity does change the value of the 
circulation—i.e., the value of I’, to place the stagnation 
point on the trailing edge—and as a result does in- 
directly change the lift. 

The increment in circulation due to the discontinui- 
ties is given by 





— B, 
ar=-)> (1m, tan <2) (30) 





(E) Net Effect of Discontinuities on Lift and Drag 


The lift and drag of an airfoil having source-sink 
discontinuities on its boundary are shown in sections 
(B) and (C) to be given by 


Lift = pUI 
Drag = pUy) ma 


Substituting for T, Eq. (29a), 





L=p Url 42RU, sin (a +e) — ) Mn tan(‘ ~ *)| 
9 
7 si (31) 
D = pUo >) m, 


The chord of this airfoil will be very nearly 4R. Thus, 
Mn e— 8B 
C, = 2r(a+6-2> — tan (“ ‘) 
n Uoc 2 
(32) 
m 
Cp =2 >» — 
” o€ 
where m,/Uoc is the flow coefficient usually denoted 
Co. 


Since, for a normal airfoil, 


C, = 2r(a + "7 
Cp = 0 (33) 


it is apparent that the increments in lift and drag due 
to the source-sink discontinuities are 





Ba — € 
= 9 7 —_———_ 
At. +2 J Co, tan ( r ) (34) 
| AC, = 2>) Co, 





The lift increment given by Eqs. (34) is plotted in Fig. 6 
vs. (8, — €)/180°. (8 — e) = O corresponds to a discon- 
tinuity at the nose, and 8 = m + e corresponds to a dis- 
continuity at the tail. Positive §’s correspond to 
points on the upper surface and negative #’s to points 
on the lower surface. The exact value of £6, corre- 
sponding to any given slot location on a specific airfoil 
is determined as in reference 3. 6, of the present 
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paper, in terms of the notation of reference 5, is given by 
B, =e =O-+6e (35) 
where 


Pay co << 
6 = sin—!V(1/2)[P + VP? + (4y/c)?] 


1 — (2x/c)? — (2y/c)? 


P = 


(x and y are measured as shown in Fig. 6 of the present 


paper) 
l 2x g- y’ 
«= -— f ¥(¢) cot dg 
2r Jo 2 


The ¢ of reference 3 will normally be small relative to 
6, so that a good approximation for 8, is obtained by 
writing Eq. (35) as 


B= e920 (36) 


The following specific conclusions can be drawn from 
these results. 

(1) A sink, or suction slot, on the upper surface of an 
airfoil increases the lift by an amount directly propor- 
tional to its strength. This increment is independent 
of angle of attack but is a function of the chordwise lo- 
cation of the discontinuity. 

(2) A sink or suction slot on the lower surface has 
just the opposite effect. 

(3) The drag increment is 2qual to the rate at which 
free-stream momentum is removed from the stream 
by the sink and is independent of angle of attack and 
position. 

(4) A source or an equivalent pressure slot has in 
every respect the opposite effect from the sink. 


(III) CoMPARISON WITH EXPERIMENT 


In order to confirm the practical validity of Eq. (34), 
a 16-012 airfoil with a single suction slot at 95 per cent c 
has been tested in the two-dimensional U.A.C. laminar 
flow tunnel. The results are summarized and com- 
pared with the theory in Fig. 8. It is apparent that the 
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effectiveness of the slot is considerably greater than 
indicated by the theory until the boundary layer has 
been removed. From here on the purely hydrodynamic 
effect of the slot becomes dominant and the rate of in- 
crease in lift is as indicated by the theory. The po- 
tential solution is thus confirmed within the limitations 
of its assumptions. 


(IV) PRacTICAL APPLICATIONS OF SOLUTION 


From a practical point of view the most interesting 
aspect of these results is the large lift increment associ- 
ated with a suction slot located at or near the trailing 
edge. Such an effect has been experimentally observed 
by Regenscheit® and is confirmed by the U.A.C. tests, 
It is to be remembered that this effect differs basically 
from that usually associated with suction slots in that 
it is a purely hydrodynamic effect and has nothing to 
do with boundary-layer removal. 

As is apparent from the solution presented in the 
previous sections, the physical explanation for this 
increase in lift lies in the effect of such a slot on the 
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A SOLUTION FOR THE LIFT AND 


trailing-edge stagnation point. The explanation is 
briefly as follows: Normally the circulation about an 
airfoil is determined by, and is just sufficient to place, 
the stagnation on the trailing edge. (See Fig. 5a.) 

When suction is applied on the upper surface and 
particularly when it is applied near the trailing edge, a 
forward increment of velocity is induced at all points 
behind the slot, including the region near the tail. As 
a result, the stagnation point occurs before the trailing 
edge is reached (see Fig. 5b). A condition analogous 
to the Kutta condition is thus established, the lower 
surface flow in attempting to flow around the sharp 
trailing edge separating. The upper surface flow ac- 
celerates in the face of the vacuum thus created and 
flows smoothly off the tail (Fig. 5c). The net effect is 
an increase in circulation, which, it can be said, is suf- 
ficient to place the stagnation point on the trailing edge. 
It is this increase in circulation which causes the in- 
crease in lift. The theory indicates that the equivalent 
of a source on the lower surface will produce the same 
effect. 

In Fig. 7 the pressure distribution on a 16-012 airfoil 
at a = 0 with a suction slot at 85 per cent is shown. 
The no-suction computation was made using the po- 
tential functions (3) and (5)* of the present report in 
The f = 0 
It is 
apparent that there is considerable (though indetermi- 
nate) external lift due to the action of the slot alone. 
But the net Lft including the internal forces as indi- 
Thus the lift 
represented by this pressure distribution is equal to‘the 
internal momentum reaction. The difference between 
the ! = 0 and I to satisfy the Kutta condition curves 
represents the net effective lift. 

It is also significant that a favorable pressure gradi- 
ent exists right up to the slot, so that such a slot would 
not only increase C;, by increasing circulation but would 
also tend to prevent separation and therefore increase 
Cumes.- 

From Fig. 7 it is apparent that, as Cg is increased, 
the stagnation point following the slot moves back to- 
ward the trailing edge. When this stagnation point 
reaches the tail and subsequently moves on under to 


combination with Theodorsen’s method. 
curve in this figure is of considerable interest. 


cated by Eq. (21) is zero for this case. 


*T is assigned the Kutta value as given by Eq. (29a). 
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the lower surface, the Kutta function of the trailing 
edge is lost, and it is not clear what will happen. In 
any event, it is clear that the analysis presented in this 
paper is no longer applicabie. Thus, although the C,; 
increment corresponding to a given Cg increases as the 
slot is moved back toward the tail, from the preceding 
comments it is apparent that the maximum Cg for 
which the analysis is applicable decreases. Thus, it 
does not necessarily follow that the maximum AC, will 
oecur with the slot as near as physically possible to the 
tail. Consideration of these points must await experi- 


mental observations. 
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Simplified Analysis of General Instability 
of Stiffened Shells in Pure Bending 


F. R. SHANLEY* 
The RAND Corporation 


SUMMARY 


Although much work!~ has been done to develop a theory for 
the failure of shells by general instability, there is at present no 
simple method by which the size of the frames may be determined 
for any given diameter, bending moment, and frame spacing. 
Such a method is needed in determining the optimum design for 
stiffened shells, to be used as a basis for weight analysis of fuse- 
lages, and other shell structures. In an extension of the work 
done for The RAND Corporation,‘ a simple coefficient has been 
determined for this purpose. Since it appears that this method 
may also be useful in design calculations, a brief description is 


presented below. 


DESCRIPTION OF THE METHOD 


bie ORIGINAL PROBLEM that caused the C.A.A. to 
inaugurate the investigations that have been under 
way for over 10 years was to determine whether the 
frames were strong enough to ensure failure of the sheet- 
stringer panels between frames—i.e., to prevent general 
instability. That is still the primary concern of the de- 
signer; we have much less need for a theory that will 
predict the ultimate bending moment of a shell in which 
the frames are not stiff enough to do this. 


For shells that are designed to fail through panel in- 
stability, it is possible to make use of the fact that a 
column loses its resistance to bending as its critical 
(buckling) load is approached. 

As the sheet-stringer elements near the most highly 
loaded point on the shell approach their critical load 
as columns, they tend to behave like hinged bars, with 
spring supports furnished by the frames.® From the 
photographs of general instability failures,* it appears 
that the proper model for a single element would be 
two bars hinged midway between two frames and at 
adjacent frames, as shown in Fig. | The two middle 
frames would be replaced by springs, for which the 
spring constant required to prevent buckling can be 
computed from the geometry of the slightly deflected 
bars. 

The equation for the required spring constant, a, 
may be written 


where P is the axial load and L is the frame spacing. 


Received May 20, 1949. 
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at Los Angeles. 


For the assumed model, A; = 3/2. For two bars, 
each of length L, supported by one spring, the con- 
stant equals 2.0. Since the value of the constant will 
not be used directly, the type of model assumed is not 
important. The only assumption that need be made 
is that the mode of failure remains the same. 

If the sheet-stringer elements are considered to be 
uniformly distributed around a shell of circular cross 
section, the maximum axial load per unit of circum- 
ference equals the bending moment divided by the en- 
closed area [see Eq. (39), reference 4]: 


mar, = 4M/rD? (2) 


where .V/ is the bending moment and D is the mean 
diameter of shell. 

If this maximum value of g is considered to apply 
over some constant fraction of the shell circumference,! 
Eq. (2) may be written 


P., = K,M/D (3) 


‘1 


This value may be substituted in Eq. (1) to obtain 
the required spring constant for the frame. The ef- 
fective spring constant for a circular frame under di- 
ametrically opposed loads may be written in the form’ 


K,(ED),/D* 4 


where (J); is the effective frame stiffness in bending. 

The value of A; may be considered to remain con- 
stant for a given distribution of the loading forces. 
(Nonlinear effects are assumed to be small and are also 
included indirectly in the final evaluation of the coef- 
ficient.) 

Eqs. (1), (3), and (4) may be combined to give an ex- 
pression for the required frame stiffness: 


(EI), = C;MD?/L (5 

If it could be shown that the coefficient C; is nearly 
constant for any given cross-sectional shape, Eq. (9 
could be used for the design of the frames. At the 
present time, the best approach seems to be to deter- 
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TABLE 1 
Test Data from N.A.C.A. Technical Note No. 909 
(Pure Bending, Aluminum-Alloy Material, E assumed = 10 X 10° Lbs. per Sq.In., P = Panel Failure, G = General Instability, P/G 


I;, In.4 
xX 105 
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. 537 
. 537 
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587 
.537 
. 5387 
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but were obtained from California Institute of Technology through Dr. 


TABLE ¢ 
Test Data from N.A.C.A. Technical Note No. 1499 
(Pure Bending, Aluminum-Alloy Material, E assumed = 10 X 10® Lbs. per Sq.In., P = Panel Failure, G = General Instability, P/G = 


Frame Spacing, J; = Moment of Inertia of Frame) 
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(See Note 1) 
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mine a value of C; from test data, such that the prob- 
ability of general instability would be reduced to an 
acceptable minimum. 

To do this, the test data from references | and 2 have 
been analyzed (ultimate bending moments for the tests 
of reference | were obtained from Dr. E. E. Sechler, 
of California Institute of Technology). (See Tables 1 
The coefficient was calculated from Eq. (5 


C; = (EI) ,.L/MD? (6 


and 2.) 


On Fig. 2, the test points have been plotted against C, 
so as to show the number of general instability and 
panel failures within each increment of C;. There is a 
fairly large range within which both types of failure 
occur. The probability of either type of failure at any 
given value of C, is more clearly indicated by the sum- 
mation curves. The top curve, for example, gives the 
number of general instability failures that occurred 
On this basis it appears that 
5 (= 1/16,000 


above a given value of Cy. 
a value of C, = 6.25 X 10 
a safe figure for the shells tested. 

The wide range of overlap does not necessarily indi- 


represents 


cate that the theory from which C; was derived is cor- 
respondingly inaccurate. In reference 2, for example, 
a wide range was given for the parameter A to describe 
the transition zone within which either type of failure 
A certain amount of such overlap will be 
In any structure designed 


may occur. 
caused by imperfections. 
on an optimum basis, the probability of a reduction in 
load capacity is increased by the fact that either of two 
types of imperfections may be expected. This also 
widens the range in which either type of failure may be 
expected to occur. 

Incidentally, it can be shown that C; can be converted 
into the form A ‘A (where A is a constant and A is the 
parameter used in reference 2) by assuming that the 
critical load per stringer is given by the Euler formula. 
One great advantage of the coefficient C; is that it does 
not contain any of the properties of the sheet-stringer 
combination. This permits the frame problem to be 
handled independently of all the usual complications of 
the panel problem, such as inelastic effects, effective 
width, interrivet buckling, etc. The coefficient C, 
brings the bending moment directly into the calculation 
of frame stiffness, while the parameter A does so in- 
directly through the effective FJ of the sheet-stringer 
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combination. It might also be possible that the critical 
conditions for general instability occur first over a fairly 
small area, in which case the present analysis would per- 
haps represent the conditions better than a theory based 
on the appearance of the structure after buckling. 

Fig. 2 shows the need for more tests in the transition 
region. It is to be hoped that future research will be 
directed toward determining how to prevent general 
instability rather than how to predict it. There is an 
urgent need, also, for tests of shells that more closely 
resemble the type-of construction actually used, with 
more efficient frames located on the inside of the shell. 

Finally, the effective -/ of various types of frames 
should be investigated, including the effects of the shell 
thickness. 
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A New High-Performance Engine Indicator 
of the Strain-Gage Type 


C. S. DRAPER* ann Y. T. LIt 
Massachusetts Institute of Technology 


ABSTRACT 


Pressure changes in the working fluid largely determine the per- 
formance of heat engines. For this reason, accurate pressure indi- 
cations are of primary importance for studies of engine operation. 
Pressure indicators have been in use for many years, but perform- 
ance defects have restricted their value even for laboratory and 
test-stand work. The operating limitations have been due to the 
fact that available pressure receivers have been sensitive not only 
to pressure but also to temperature, to mechanical strains due to 
mounting effects, and to mechanical vibrations in the walls of the 
pressure chamber. It is the purpose of this paper to describe the 
features of a new pressure receiver using a high-flexibility cate- 
nary-type diaphragm and a temperature-compensating wire- 
wound strain gage to give satisfactory pressure-measuring char- 
acteristics with greatly reduced response to temperature and ex- 
traneous mechanical effects. 


(I) PERFORMANCE REQUIREMENTS FOR ENGINE 
PRESSURE INDICATORS 


ee ee INDICATORS receive pressure vari- 

ations in the working fluid and produce indications 
or records of these fluctuations. Records of two general 
types are in common use: (1) average pressure records 
and (2) continuous pressure records. 

Average pressure records are possible only from equip- 
ment that operates with repetitive cycles, such as re- 
ciprocating piston engines. Continuous pressure rec- 
ords may be taken from all kinds of engines and tur- 
bine-type power plants and from research chambers us- 
ing nonrepetitive explosions. 

To be generally satisfactory, continuous pressure in- 
dicators!-*5, °% 31 must be able to withstand engine op- 
erating conditions for reasonable periods and give sat- 
isfactorily accurate records for both low- and high- 
frequency pressure variations. The low-frequency 
variations to be measured are directly associated with 
the power generation phase of engine operation; the 
high-frequency variations involved may be associated 
with detonation, shock waves inside a closed chamber, 
or air pulsations within an induction system. In addi- 
tion to these primary operating characteristics, the 
pressure datum line and the calibration of a satisfactory 
indicator must be unaffected by environmental effects 
such as temperature and extraneous strains. Unless 
this condition is fulfilled, pressure measurements of 
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high accuracy are difficult to realize, because a calibra- 
tion based on a constant environmental temperature 
will not be valid for an actual indicator installation sub- 
jected to the fluctuating temperatures accompanying 
power-plant operation. 

In addition to the calibration-changing effects of tem- 
perature and installation strains, pressure indicators 
are subject to dynamic errors due to inertia reaction 
forces acting on mechanical parts. These forces pro- 
duce distortions in the records of rapidly changing pres- 
sures and, in addition, cause an indicator to respond to 
vibrations of the wall to which it is attached. The re- 
duction of inertia reaction effects to tolerable levels is 
one of the primary objectives of pressure-indicator de- 
sign. 

Averaging indicators*®-** are designed to overcome 
the inaccuracy-producing effects usually present in con- 
tinuous pressure indicators by removing the calibrated 
parts of the system from the influence of engine tem- 
perature and vibration. This is achieved by designing 
the system so that only a pressure equality sensitive unit 
is exposed to the working fluid pressure. In the indi- 
cator, a slowly varying pressure is applied to a cali- 
brated pressure recording unit, and each engine cycle 
places a single point on paper carried by a rotating 
drum connected to the engine shaft, when the working 
fluid pressure becomes equal to the external pressure. 
This system satisfactorily avoids inaccuracies due to 
temperature and dynamic effects but can be applied 
only to engines giving repetitive cycles and, in any case, 
is unable to record high-frequency pressure oscillations. 

The discussion that follows describes a continuous 
pressure indicator that is substantially free from the 
common defects of previously existing indicators. 
Design features are described, and the performance 
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characteristics realized are reported from data taken 
in laboratory and engine tests. 


(Il) GENERAL FEATURES OF ENGINE INDICATORS 


Piston and Spring Engine Indicator* 

Engine pressure records were first made by means of 
the piston and spring arrangement introduced by James 
Watt and commonly called the steam engine indicator. 
A schematic diagram of this arrangement is shown in 
Fig. 1 and a corresponding functional diagram is given 
in Fig. 2. In this type of indicator, the part of the pis- 
ton exposed to the pressure to be measured supplies the 
pressure seal, which receives pressure as the actuating 
input and applies a corresponding force to the piston, 
which acts as the force summing member. A coil spring 
placed between the body of the indicator and the piston 
is the elastic restraint, which reacts on the force sum- 
ming member with a force proportional to the linear 
displacement of this member with respect to the cylin- 
der. A mechanical linkage serves as the coupling system 
by transferring the pressure receiver output deflection 
to the indicating system, which consists of the recording 
index and a rotating drum driven by the engine shaft. 

The usefulness of the steam engine indicator is lim- 
ited by the relatively large mass of its moving parts and 
the amount of gas flow required to change the pressure 
acting on the piston. These characteristics make it 
impossible for conventional piston and spring indicators 
to give satisfactory records of the high-frequency oscil- 
lations associated with detonation, or even the rapid 
pressure fluctuations of high-speed engine operations. 


Averaging Pressure Engine Indicator™®~” 


The schematic diagram of Fig. 3 illustrates the prin- 
ciple used by designers of balanced pressure diaphragm 


FUNCTIONAL DIAGRAM FOR PRESSURE INDICATOR USING A PRESSURE RECEIVER WITH MECHANICAL DISPLACEMENT OUTPUT 


indicators to retain the simplicity and reliability of 
the piston and spring system for converting pressure 
levels into indicator record deflections and, at the same 
time, to remove calibrated parts from the severe tem- 
perature and vibration conditions of engine operations. 
As shown in this figure, an external source of high pres- 
sure is applied to the piston and to the outer side of a 
diaphragm exposed to the engine pressure on its inner 
surface. A contact that is effective in the primary cir- 
cuit of an induction coil causes a spark to pass from the 
recording index electrode through recording paper when 
the pressure moves the diaphragm against or away from 
the contact to give a make or break of the circuit.* 
The recording paper is carried by a rotating drum 
driven from the engine shaft. An indicator card is 
made by slowly varying the external pressure until the 

* In practice, an electronic circuit is often used between the 
pressure equality indicating diaphragm and the induction coil 
With this arrangement, the current broken by the contacts ts 
greatly reduced from the current that would be required if the full 
primary current passed through the contacts. 
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FIG.3 SCHEMATIC DIAGRAM 0 
AVERAGING PRESSURE INDICATOR 
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single points, or pair of points, made on the record for 
each power cycle of the engine combine to produce a 
complete average pressure crank angle record. 

In addition to the system of Fig. 3, average-pressure 
records may be made up from a piston, spring, and drum 
arrangement like that shown in the figure but including 
a sampling valve*® to replace the balanced pressure dia- 
phragm. 
pressure vessel wall, is arranged to open for a short time 


The sampling valve, which is mounted in the 


at a controllable crank position angle corresponding to 
the recording drum position. This causes the pressure 
acting on the indicating piston and spring to build up to 
a level substantially equal to the engine-cycle pressure 
for the crank angle setting used. By making settings 
for the proper number of crank angles, a complete indi- 
cated pressure cycle record may be taken. 

It is obvious from the description of the preceding 
paragraph that averaging indicators are not suitable for 
making pressure records of rapidly fluctuating pressures 
that are not exactly repeated a large number of times. 
Because of this defect, averaging indicators cannot be 
used for recording nonrepetitive pressure changes or 
high-frequency oscillations. 


Flat-Diaphragm Type Engine Indicator: ® * ®, 5-11, 13, 31 


The need for continuous pressure records in all funda- 
mental investigations of engine performance has pro- 
duced many attempts to construct equipment without 
the operating deficiencies of the piston and spring in- 
dicator. In this direction, one of the fundamental 
modifications to the system of Fig. 1 has been to reduce 
the ratio of the mass to the elastic restraint of the mov- 
ing parts of the pressure receiver in order to improve its 
response to rapid pressure changes. In many designs, 
this has been achieved by replacing the piston and 
spring system with a flat diaphragm, as indicated in the 
schematic diagram of Fig. 4. In this diagram, the flat 
diaphragm combines within itself all the functions of 
the pressure receiver of Fig. 2 in the sense that it acts as 
the pressure seal, the elastic restraint, and the force 
summing member. At the same time, it determines 
the equivalent mass of the moving system and supplies 
the pressure receiver output as the effective displace- 
ment of its central area 

The output displacements of a flat diaphragm pres- 
sure receiver are generally so small that the coupling 
and indicating systems used must be based on either 
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NGINE INDICATOR 595 
optical'* *! or electrical’ principles. Because of their 
fragility and bulk, optical systems are more useful for 
special laboratory investigations than for general en- 
gine work. When, as in Fig. 4, the diaphragm output 
displacement is used to vary an air gap, the electrical 
coupling system can be designed to operate on the re- 
sulting change in either the electromagnetic reluctance! 
or the electrostatic capacity.” '' In an indicator of this 
kind, an electronic amplifier completes the coupling 
system and supplies electrical signals to some type of re- 
cording oscillograph. 

At the present time, amplifiers and oscillographs are 
so well developed that the electrical components re-~ 
quired for high-performance coupling and indicating 
systems are easily available. There are many designs 
possible for systems of this type which will operate 
without imposing functional limitations of any kind on 
the pressure-indicating equipment in which they are 
used, for which reason the performance possible from a 
pressure indicator is usually determined mainly by the 
characteristics of the pressure receiver. 
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(III) PERFORMANCE LIMITING CHARACTERISTICS OF 
FLAT-DIAPHRAGM PRESSURE RECEIVERS 


General Characteristics of Flat-Diaphragm Pressure 


Receivers'~* 

Engine indicators intended for use with high-speed 
engines have usually been designed with pressure re- 
ceivers that are either flat diaphragms or have the flat- 
diaphragm shape nearly enough for the elastic theory 
associated with this shape to be valid within reasonable 
limits. For this reason, it is possible to discuss the gen- 
eral problems of pressure-receiver design in terms of the 
. elastic theory of flat diaphragms.* * 

The schematic diagram and the corresponding func- 
tional diagram, shown in Fig. 5, represent the essential 
operating characteristics of a flat-diaphragm pressure 
receiver under static conditions. By definition, static 
conditions exist when both the essential input and the 
corresponding output are constant. The static char- 
acteristics can be expressed in terms of sensitivities, 
where a particular sensitivity is determined as the ratio 
of a small output change to the corresponding small in- 
put change when all the other essential quantities are 
held constant. 
sensitivity of a flat diaphragm pressure receiver is equal 
to the ratio of the displacement of the diaphragm center 
to the pressure applied to the diaphragm. 


For example, the pressure - displacement 


Metallic diaphragms applied as pressure receivers are 
affected by four input quantities: 

(1) Pressure applied to the outer surface of the dia- 
phragm. 

(2) Linear acceleration of the diaphragm support in 
the direction at right angles to the diaphragm surface. 

(3) Temperature gradient in the material of the dia- 
phragm between its inner and outer surfaces. 

(4) Average temperature of the diaphragm material. 

Pressure, the first of these quantities, is the desired 
actuating input for a pressure receiver. Linear accelera- 
tion of the diaphragm support and temperature gradi- 
ent in the diaphragm material cause undesired deflec- 
tions that are not distinguishable from pressure effects. 
Changes in the average temperature of the diaphragm 
do not of themselves cause output deflections but act as 
a modifying input by affecting the elastic constants 
(modulus of elasticity and Poisson's ratio) of the dia- 
phragm material and so altering the sensitivities of the 
diaphragm for other inputs.* 4 

In addition to its static performance, as measured by 
the sensitivities of Fig. 5, the dynamic characteristics 
of a pressure receiver are important. By definition,’ 4 
dynamic conditions exist when the relationship be- 
tween an output quantity and the corresponding input 
quantity depends upon the rates at which these quanti- 
ties are changing. For an operating component like 
the flat diaphragm, which has relatively low damping, 
the essential factor in determining the dynamic char- 
acteristics of the component is its undamped natural 
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DIAPHRAGM MATERIAL: STEEL 
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OF FLAT DIAPHRAGM 


frequency. Fig. 6 gives a plot of the undamped natural 
frequency for a '/2-in. diameter flat steel diaphragm asa 
function of the diaphragm thickness.* * It appears 
from this plot that the undamped natural frequency of 
this flat steel diaphragm is a linear function of the dia- 
phragm thickness, with a thickness of 0.064-in. giving a 
natural frequency of about 100,000 cycles per sec. In 
general, for a given frequency of the input pressure, the 
undesirable effects of dynamic conditions are reduced 
as the natural frequency of the diaphragm is increased. 
The second curve of Fig. 6 shows the considerable cost, 
in terms of reduced pressure - deflection sensitivity, of 
increasing the diaphragm thickness to improve the 
dynamic response characteristics. 

In practice, the coupling system - indicating system 
combination of a pressure indicator receives the output 
deflection of the pressure receiver and produces a de- 
flection of the recording index. The input deflection is 
received without consideration of the extent to which 
the effects of linear acceleration and temperature gra- 
dient may have entered into the resultant. Because the 
indicating system record is interpreted in terms of pres- 
sure, it is important to have information available for 
estimating the effects of linear acceleration and temper- 
ature gradient in terms of the pressure acting on the 
pressure receiver. This information is supplied by the 
sensitivity ratios defined in Fig. 5. The acceleration- 
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pressure sensitivity ratio is the ratio of the acceleration - 
deflection sensitivity to the pressure - deflection sensitiv- 
ity. It has the dimensions of pressure divided by ac- 
celeration and represents the amount of pressure that 
would be required to produce the same diaphragm out- 
put deflection as a unit acceleration of the support. In 
asimilar way, the pressure - temperature gradient sensi- 
tivity ratio represents the pressure that would be re- 
quired to produce the same output deflection as a unit 
temperature gradient through the diaphragm. 


Effect of Linear Acceleration of the Support on a Flat- 

Diaphragm Pressure Indicator 

The curves of Fig. 7 show the way that the accelera- 
tion - pressure sensitivity ratio varies as a function of 
the ratio of the diaphragm support forcing frequency to 
the undamped natural frequency of the diaphragm for 
various diaphragm thicknesses. In an ideal pressure 
receiver, acceleration of the support would have no 
eflect on the output displacement; the corresponding 
curve in Fig. 7 would be a straight line along the fre- 
quency ratio axis. The actual curves of Fig. 7 show that 
the sensitivity ratio value of a flat-diaphragm pres- 
sure receiver increases as the diaphragm thickness (and, 
consequently, the diaphragm natural frequency) is 
made larger and changes toward higher sensitivity ratio 
values with increasing slope as the frequency ratio in- 
creases. 
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FiG.6. TEMPERATURE DIFFERENCE INPUT 
DIAPHRAGM TYPE PRESSURE RECEIVER AS A FUNCT 


A requirement of the usual pressure-indicator design 
problem is that the undamped natural frequencyeof the 
pressure receiver must be as great as, or greater than, 
some lower limit depending upon a specified inaccuracy 
tolerance for the highest pressure oscillation frequency 
to be recorded. The curve family of Fig. 7 shows that a 
flat-diaphragm pressure receiver has the characteristic 
of a greater response to acceleration, in terms of pres- 
sure, as the diaphragm thickness is increased. It fol- 
lows that, once a flat diaphragm has been made suffi- 
ciently thick to withstand the rigors of operating condi- 
tions and to have the minimum required undamped nat- 
ural frequency, the undesirable effects of support ac- 
celeration will always be made worse with respect to the 
effective pressure-controlled output by increasing the 
diaphragm thickness. This means that the flat-dia- 
phragm pressure receiver has a limitation so far as con- 
cerns the accurate recording of pressure fluctuations 
that have the same frequency as the imposed support 
vibration. When detonation occurs in an engine cylin- 
der, the walls of the cylinder usually have vibratory ac- 
celeration frequencies equal to the existing detonation 
frequencies, so that records from a pressure indicator, 
using the flat diaphragm as,its pressure receiver, will be 
subject to inaccuracies that are difficult to eliminate 
by compensation and practically impossible to remove 
effectively by corrections to the actual records. 
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Effect of Temperature Gradient in the Material of a Flat- 

Diaphragm Pressure Receiver 

The theoretical curve of Fig. 8 is practically equiva- 
lent t8 a curve of temperature difference - pressure sen- 
sitivity ratio for a flat - diaphragm pressure receiver as a 
function of the diaphragm thickness. For convenience 
in presentation, the reference quantity is taken as the 
temperature difference between the outer and inner sur- 
faces of the diaphragm, instead of the temperature gra- 
dient itself. The curve shows that the temperature dif- 
ference - pressure sensitivity ratio is generally negative. 
This means that a positive temperature difference (the 
outer-surface temperature being taken as the reference 
temperature) produces the same effect as a reduction of 
the input pressure. Physically, the action is caused by 
the greater thermal expansion of diaphragm layers near 
the hotter surface—that is, when the outer surface is 
hotter than the inner surface, the diaphragm tends to 
bow outward and produces the same effect as a reduc- 
tion in the applied pressure. This effect is considerable, 
as illustrated by the fact that, with a diaphragm of 
0.030-in. thickness, a 1°F. temperature difference pro- 
duces the same effect as a pressure change of about 10 
Ibs. per sq.in. In practice, the temperature gradient in 
a pressure receiver exposed to the operating pressures of 
an engine cylinder will change-considerably as engine 
conditions are changed. As a consequence, the effec- 
tive zero reference level on the record from a flat - dia- 
phragm type pressure receiver varies in a way that usu- 
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ally introduces inaccuracies of considerable size in the 
corresponding engine pressure records. 

In addition to the temperature gradient effect due to 
a change in engine conditions, there is a more serioys 
effect caused by the change in temperature that occurs 
within a single cycle. Because the cyclic temperature 
effect and the pressure input are of the same fundamen. 
tal frequency, there is no way to distinguish one from 
the other in the output signal of the pressure receiver. 


Effect of Average Temperature on the Performance of Flat- 

Diaphragm Pressure Receivers 

Fig. 9 gives a plot of the pressure - displacement sensi- 
tivity performance ratio for a flat diaphragm pressure re- 
ceiver as a function of the average diaphragm tempera- 
ture. The performance ratio plot shows that the pres- 
sure - displacement sensitivity increases as the average 
diaphragm temperature becomes greater. This change 
corresponds to the physical action that causes a weaken- 
ing of the diaphragm material as its average temperature 
increases. The effect represented in Fig. 9 does not 
cause a change in the zero reference level of a pressure 
receiver, but, on the other hand, it introduces an inac- 
curacy in the interpretation of pressure records if the 
average diaphragm temperature is either at a constant 
unknown level or fluctuates in an unknown way. In en- 
gine installations, the average diaphragm temperature 
has been known to vary several hundred degrees 
Fahrenheit;* this effect can introduce inaccuracies of 
several per cent in engine-pressure records. 
Pressure 


Performance Defects of a Flat-Diaphragm 


Receiver'-* 

Specifications for the operating components of any 
measuring system can be written only after acceptable 
inaccuracy tolerances have been established for indica- 
tions from the system. In the case of pressure indicators, 
it appears reasonable to establish an inaccuracy toler- 
ance equal to 1 per cent of the designed maximum pres- 
sure indication. Even though no inaccuracy compo 
nent is allowed for the coupling and indicating systems, 
the preceding discussion of the flat diaphragm shows 
that pressure receivers of this type have performance 
characteristics that may produce inaccuracy compo 
nents larger than the specified tolerance. For example 
when a steel diaphragm of !/»-in. diameter and 0.050- 
thickness is assumed for the pressure receiver of an indi 


TABLE 1 
Tolerance on Environmental Effects Equivalent to Pressure E 
fects of 10 Lbs. per Sq. In. on Flat Steel Diaphragm 0.050 It 









Thick and 0.500 In. in Diameter 
Support Vibra- 
tion Ampli- 
tude at 10,000 
Cycles per Sec.; 
(Based on Curve 
of Fig. 7) (In.) 


Temperature Dif- 
ference Through 
Diaphragm; 
(Based on Curve 
of Fig. 8) (°F.) 


erage Tempef 
ture When 4p 
plied Pressure! 
500 Lbs. per % 
In. (Based & 
Curve of Fig.’ 
(°F.) 


- 
25 


1.04 X 10-* 





Difference in Av 






















an 
ab 


ati 


rec 


tu 
en 
en 
ere 
dif 
ina 
alle 


ma 
ent 
in { 
rece 
cee 
tim 





e in the 


t due to 
> serious 
t occurs 
perature 
idamen- 
ne from 
eiver, 


2 of Flat- 


nt sensi- 
ssure re- 
empera- 
he pres- 
average 
s change 
weaken- 
perature 
loes not 
pressure 
an inac- 
ds if the 
constant 
y. In en- 
perature 

degrees 
racies of 


Pressure 


s of any 
ceptable 
yr indica- 
dicators, 
cy toler- 
um pres- 
- compo- 
systems, 
mi shows 
‘ormanice 
- Compo 
example, 
().050-in. 
f an indi 


-essure El 
. 0.050 In 


re in AY 
Temper 
When Ap 
Pressure» 
bs. per % 
Based @ 
of Fig. ! 
F.) 

5 
















A NEW 


HIGH-PERFORMANCE 





ENGINE INDICATOR 




































































eoocococor rrr rrr eee eee COUPLING SYSTEM ~~~ mem mmm me ae ee ee ee > 
' ' 
(cs) 
l l COUPLING 
| , System 
PRESSURE RECEIVER ELECTROMAGNETIC CHANGE ouTPuT 
ouTPuT | oR + IN RELUCTANCE ELECTRICAL 1 ont 
DISPLACEMENT | EFFECTIVE CHANGE IN GAP ELECTROSTATIC OR CAPACITANCE TRANSDUCER | 
Rei I Xo=Xor + Xism CIRCUIT A(rel) or A(cap) (et) 1 Gs 
(emc, on esc) 
ba 
z 
1 wt } 1 
| wy ! 
ave | 
| auto | 
ad2z ! 
| 7) rs = | 
| ous ! | 
‘ a6 | 
AVERAGE | ag | 
TEMPERATURE | | . I 
Tisr | SHELL STRUCTURAL x i 1 
> | DEFORMATION TOTAL GAP DISTANCE ai | 
NSTALLATION ee 
srmess | | (Sh) Xion gd=X, +(9d), 
OF SHELL | 
! 
(is).  ! i 
| \ 
! I 
Ihc a a eae shes Sins “ari ied putin mashed aes as ec mee na tis die ies sas ons ea i a ees db ite hae dien J 


FIG. 10. FUNCTIONAL DIAGRAM OF DISPLACEMENT RECEIVER FOR ELECTRICAL COUPLING SYSTEM. 


cator designed to measure a maximum pressure of 1,000 
lbs. per sq. in., Table 1 gives the allowable levels of vi- 
bratory acceleration, temperature difference across the 
diaphragm, and average temperature change to produce 
an output deflection change equivalent to a pressure 
variation of 10 Ibs. per sq.in. (i.e., 1 per cent of 1,000 
Ibs. per sq.in.). 

Available data on cylinder head vibration‘ show that 
amplitudes of 0.0002-0.0005 in. occur at frequencies 
above 10,000 cycles per sec. This means that the flat 
steel diaphragm does not have a sufficiently low acceler- 
ation - pressure sensitivity ratio to serve as the pressure 
receiver of a satisfactory engine indicator. 

Accurate information is not available on the tempera- 
ture differences through steel diaphragms exposed to 
engine pressure, but it is reasonable to expect differ- 
ences up to 100°F. under operation conditions.’ Ref- 
erence to Table 1 shows that the effect of a temperature 
difference through the diaphragm can cause a pressure 
inaccuracy that exceeds by a considerable margin, the 
allowable limit of 10 Ibs. per sq.in. 

The average temperature of a flat steel diaphragm 
may change from about 70° to about 300°F. with differ- 
ent operating conditions. The corresponding variation 
in the output deflection of a flat diaphragm pressure 
receiver for an applied pressure of 500 Ibs. per sq.in. ex- 
ceeds the specified tolerance limit of Table 1 by ten 
times. 


Performance Limitin g Characteristics of a Flat Diaphragm 
Pressure Receiver Electrical Displacement Detecting 
Coupling System 
The combination of a flat diaphragm pressure re- 

ceiver and an electrical deflection detecting coupling 

system, as shown in Fig. 4, is most commonly adopted 
in indicator design because of the simplicity and rug- 


gedness of the system and the good performance in low- 
frequency applications.® 5 


However, performance 





limiting characteristics, in addition to those already de- 
scribed for the flat diaphragm as a pressure receiver, are 
introduced by the use of the electrical displacement de- 
tecting coupling system in which the flat diaphragm 
serves also as an element of the electrical system. 

An electrical displacement detecting coupling system 
of either the electromagnetic or the electrostatic type re- 
ceives as the input signal the diaphragm displacement, 
which controls the reluctance or capacitance across an 
air gap between the diaphragm and a “‘fixed’”’ pole or 
electrode, and thereby produces an electrical signal that 
is amplified as the output signal of the coupling system. 
As shown in Fig. 4, the flat diaphragm pressure re- 
ceiver and the fixed pole or electrode are linked together 
by the outer shell of the indicator, and a diaphragm de- 
flection is realized as an effective change in gap only 
through the bridging of the structure of the shell. The 
functional diagram of the essential operating character- 
istics of the deflection detecting coupling system is 
shown in Fig. 10. Under static conditions, the effec- 
tive change in gap between the two electrodes (one of 
which is the flat diaphragm) are affected by three input 
quantities: (1) deflection of the diaphragm as a pres- 
sure receiver; (2) deformation of the shell caused by 
temperature change; and (3) deformation of the shell 
caused by externally applied stresses due to mounting 
effects. 

The second and third of the foregoing inputs are un- 
desired; their existence produces a shift in the datum 
line of the pressure record. In actual test of a typical 
design, the zero pressure datum line shift of the pressure 
record due to each of the two effects may, under ordi- 
nary operating conditions, exceed the record reading cor- 
responding to the maximum allowable pressure of the 
diaphragm. 

In Fig. 10, the effective change.in gap is shown as the 
controlling input of an electromagnetic or an electro- 
static circuit, the output of which is a change in induct- 
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ance or capacitance. Under ideal conditions, the gap 
change input - reluctance (or capacitance) change output 
sensitivity of the circuit should be a constant. However, 
as an inherent characteristic of the circuit, the gap 
change input - reluctance (or capacitance) change output 
sensitivity is modified by the absolute value of the total 
gap. Consequently, the diaphragm deflection will not 
only provide an input signal to the electromagnetic (or 
electrostatic) circuit but will also impose a modification 
upon the sensitivity of the circuit through the change in 
the total gap. A nonlinear sensitivity is therefore an in- 
evitable result. 

By a similar reasoning, the temperature and installa- 
tion-strain effects applied to the shell of the indicator 
will not only produce a shift of the zero pressure datum 
line in the record but will also cause the gap change in- 
put - reluctance change output sensitivity to undergo a sub- 
stantial alteration. Actual tests of a flat diaphragm - 
electrical displacement detecting type indicator show a 
deviation from unity of as much as 10 per cent in the 
sensitivity ratio because of temperature and installa- 
tion strain effects. 

The foregoing discussion shows that the flat dia- 
phragm does not have characteristics that permit it to 
be satisfactory as the receiver of a pressure indicator. 
When flat-diaphragm performance is used as the stand- 
ard of comparison, it appears that a suitable pressure 
receiver must have a reduced sensitivity to acceleration 
and temperature, in addition to having a pressure - dis- 
placement sensitivity that is less affected by average 
temperature. The essential features of a pressure re- 
ceiver substantially free from the performance defects 
of the flat diaphragm are described in the following sec- 
tion. 


(IV) THE CATENARY DIAPHRAGM-STRAIN 
GENERATING TUBE COMBINATION PERFORMANCE 
CHARACTERISTICS FOR PRESSURE RECEIVERS 


Desired Features of a Satisfactory Pressure Receiver 


The problem of producing a pressure receiver with 
sufficient mechanical strength and temperature re- 
sistance to reliably withstand engine operating condi- 
tions can be solved by separating the pressure seal from 
the member that is designed to provide elastic reaction 
to balance the force output of the pressure seal. This 
feature makes it possible to design the pressure seal to be 
so weak as an elastic member and so small in effective 
mass that temperature and acceleration effects do not 
appreciably change the resultant output of the pressure 
receiver. The member required to complement the elas- 
tically weak pressure seal must be strong from the 
standpoint of elastic deformations and must have a 
high natural frequency even when loaded by the mass of 
the pressure seal. In addition, it is desirable for the 
quantity chosen as the pressure receiver output to be 
generated in a way that will make it possible to incor- 
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FIG. 11. PICTORIAL SCHEMATIC DIAGRAM OF _CATENARY DjiA. 
PHRAGM - STRAIN GENERATING TUBE COMBINATION PRESSURE &€. 
CEIVER WITH WIRE STRAIN GAGE WINDING AS RECEIVER FOR ELECTRICAL 
COUPLING SYSTEM. 


porate acceleration- and temperature - compensating | 


features in the design for the system. 


The desired condition of a weak pressure seal is par- 
tially fulfilled in the design of the carbon-resistor-type 
indicator® or a piezoelastic-type indicator.** In each, a 
relatively thin flat diaphragm is used as the pressure 
seal, and the major elastic restraint of the indicator is 
supplied by a member not in direct contact with the 
hot gas. Indicators of this type may have better per- 
formance than those which have unsupported flat dia- 


phragms. However, with the pressure seal diaphragm 





only partially supported, temperature and mounting | 
strain still have important effects on the performance of | 


the indicator. 


Essential Features of a Catenary-Type Diaphragm 


Fig. 11 is a pictorial schematic diagram showing the 
essential features of a pressure receiver made up of a 
catenary-type diaphragm, which bears lightly against a 
force summing member that has the form of a rounded 
circular ring on the unsupported end of a thin tube rig- 
idly attached to the case of the unit. The catenary dia 
phragm is designed to be strong enough, as a pressure 
seal, to support the input pressure and, at the same time, 
to be so weak flexurally that any deflection occurring it 
practice will not cause an elastic force of appreciable 
magnitude to be applied to the force summing met 
ber. Under static conditions, the pressure seal trans 
forms the input pressure into a force, applied to the 
force summing member, which is substantially balanced 
by a reaction force generated by strains in the strait 
generating tube. 


Structurally, a three-dimensional catenary diaphragm 


and strain generating tube system is closely analogous to! 


(two-dimensional) suspension bridge with catenary: 


shaped cables. A similar analogy links a flat diaphragm ff 
to a solid beam. A suspension bridge is known to bf 
structurally more efficient than a solid beam bridge F 
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similarly, a catenary diaphragm - strain generating tube 
combination is more efficient in its structural design 
thana flat diaphragm. To its credit, a smaller effective 
mass is added to the catenary diaphragm - strain gener- 
ating tube system than is added in a flat diaphragm of 
equal elastic properties. As a consequence, the inertia 
force experienced by the moving parts of the pressure 
receiver during a dynamic response is less with a cate- 
nary diaphragm -strain generating tubesystem than with 
a flat diaphragm. 

Fig. 12 gives a line schematic diagram, a functional 
diagram, and a summary of static performance charac- 
teristics for the catenary diaphragm - strain generating 
tube pressure receiver of Fig. 11. The essential output 
of the pressure receiver is the resultant surface strain of 
the strain generating tube. Sensitivities for which pres- 
sure, acceleration, and temperature gradient through 
the catenary diaphragm are the respective actuating in- 
puts and resultant surface strain of the strain generating 
tube is the output are summarized in Fig. 12. 
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Improvement in Acceleration - Pressure Sensitivity Ratio 


The improvement in reduced sensitivity to the sup- 
port acceleration, which is realized by using the cate- 
nary diaphragm -strain generating tube combination in- 
stead of the flat diaphragm to carry out the functions of 
a pressure receiver, is shown by the curve family of Fig. 
13. These curves, which are drawn for various lengths 
of the strain generating tube, are similar to the curves of 
Fig. 7 for the flat-diaphragm pressure receiver. Typical 
flat-diaphragm curves from Fig. 7 are included in Fig. 13 
for comparison purposes. It appears that the accelera- 
tion - pressure sensitivity ratio for a catenary dia- 
phragm - strain generating tube combination is consider- 
ably smaller than the acceleration - pressure sensitivity 
ratio for the flat diaphragm. For a given thickness of a 
flat diaphragm and a given length of a strain generating 
tube to give comparable natural frequencies (for exam- 
ple, 0.050-in. thickness and 0.500-in. tube length), the 
relative reduction in acceleration response becomes 
greater as the frequency ratio (i.e., the forced frequency 
of the support vibration) increases. A second advan- 
tage of the catenary diaphragm - strain generating tube 
combination is that the acceleration - pressure sensitiv- 
ity ratio for a given forcing frequency decreases as the 
natural frequency of the combination is increased. The 
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FIG 14. FUNCTIONAL DIAGRAM FOR PRESSURE INDICATOR USING CATENARY 


curves of Fig. 13 show that, for pressure receivers with 
undamped natural frequencies between 75,000 and 80,- 
000 cycles per sec. that are,subjected to forcing fre- 
quencies near 10,000 cycles per sec. (this is in the fre- 
quency range of cylinder wall vibration of practical en- 
gines), the equal support acceleration amplitudes will 
produce from the catenary diaphragm unit an output 
strain component of only about one-fifth of the pressure 
equivalent output that would be caused by the same vi- 
bratory acceleration applied to a flat-diaphragm unit. 


Improvement in Effect of Temperature on Pressure-Strain 

ensitivity 

From the standpoint of theory, temperature gradi- 
ents in a sufficiently flexible catenary diaphragm will 
cause forces of negligible magnitudes on the force sum- 
ming member in comparison with the force output of 
the pressure seal. Any catenary diaphragm - strain gen- 
erating tube combination for which this is true will 
show no change in its pressure - strain sensitivity as the 
average diaphragm temperature is changed, because a 
variation in an already negligible quantity (the flexural 
stiffness of the diaphragm) cannot have an appreciable 
effect on the output. Similarly, a deformation in the 
outer shell of the indicator produced by a change in its 
temperature or by a change in the mounting stress will 
have a negligible effect upon the axial loading of the 
strain-generating tube, because the catenary diaphragm 
is substantially flexible. For these reasons, a properly 
designed catenary diaphragm -strain generating tube 
combination pressure receiver is theoretically free from 
temperature and installation strain effects. The reali- 
zation of this ideal performance in an actual pressure 
receiver depends upon the design of proper temperature 
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ESSENTIAL ELEMENTS OF PRESSURE RECEIVER 
TUBE COMBINATION AS THE PRESSURE RECEIVER 


compensation into the strain-receiving component of 
the coupling system. The strain-receiving means ac- 
tually used and its performance as a component of a | 
complete pressure indicating system are discussed in the 
remaining sections. 





(V) CoMPLETE PRESSURE INDICATOR USING 
CATENARY DIAPHRAGM - STRAIN GENERATING TUBE 
COMBINATION PRESSURE RECEIVER 


Functional Description of Complete Pressure Indicator 


Fig. 14 gives the functional diagram for a complete 
pressure indicator that is based on the catenary dia- 
phragm pressure receiver. The features of the pressure 
receiver have already been discussed. In effect, this op- 
erating component receives the input pressure and pro- 
duces surface strain of the strain generating tube as its 
output. The strain receiver for the coupling system isa 
temperature-compensating wire strain gage wind 
ing.7"3?-34 The output of this strain gage is received by 
the coupling system amplifier, which produces the cou- 
pling system signalasitsoutput. The indicating system 
receives the coupling system output signal as its input 
and produces the indicating system index deflection as 
its output. The pressure scale, which may be attached 
to the indicating system or may be applied separately 
to individual records, has the function of showing the 
correspondence between indicating index displacemet! 
and measured pressure levels. 


Coupling System 


The function of the coupling system is to receive the 
strain signal output of the pressure receiver and to pro 
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duce as its output an electrical signal suitable for use as 
an input by the indicating system. To be satisfactory, 
the coupling system output must be a linear function of 
the strain input without being affected by pressure re- 
ceiver temperature or accelerations. This result is ac- 
complished by using an electrical winding, with two 
arms of equal resistance, attached in a lattice form to 
the outside of the strain generating tube by a special 
temperature-resisting cement. As shown in the dia- 
gram of Fig. 15, one section of the twin-arm winding is 
applied with the wire cemented parallel to the longitu- 
dinal axis of the strain generating tube and running 
from end to end of the tube in loops equally spaced 
around its circumference. The other winding section 
is applied as a double spiral around the circumfer- 
ence of the tube in a bifilar arrangement to minimize 
the inductance effect. 

In an electrical circuit these two windings are con- 
nected in series, with a known voltage applied across 
them. The potential change at the common lead of the 
two windings is used as the output signal of the strain 
receiver. When an axial load is applied to the strain 
generating tube, the longitudinal winding is subjected 
toa reduction in electrical resistance because of com- 
pressive strain; at the same time, the cross-wound coil 
experiences an increase in resistance because of tensile 
strain. These effects are additive and cause a change in 
potential at the center tap of the series winding which 
establishes an output signal corresponding to the input 
pressure applied to the diaphragm. 

Because both windings of the strain receiver are inti- 
mately bound to the strain generating tube, they closely 
follow temperature changes in the tube and receive sim- 
ilar resistance changes due to a change in the tempera- 
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OF CATENARY DIAPHRAGM 


PRESSURE RECEIVER 


ture of the tube and wire. As a result, no potential 
change will occur at the output terminal (i.e., the com- 
mon lead of the two strain receiver windings) because of 
temperature effect. 

In addition to this, the pressure input-resultant strain 
output sensitivity holds a constant value, despite changes 
in the temperature of the strain generating tube, because 
of the compensation effect between the change in the 
modulus of elasticity and the change in Poisson's ratio 
of the strain generating tube. With only the longitudi- 
nal strain as the output signal, the pressure-strain sensi- 
tivity value at higher temperature would be larger than 
at lower temperature because of a change in the modulus 
of elasticity. However, the circumferential strain ex- 
periences an opposite effect because of the change in 
Poisson's ratio. The two effects just about balance 
each other, and, as a result, the sensitivity performance 
ratio is nearly independent of temperature, as shown in 
Fig. 16. 

Fig. 17 shows the internal electrical connections of 
the coupling system that receives resistance changes 
from the strain receiver and produces electrical signals 
suitable for operating the indicating system. The elec- 
trical system consists of a voltage supply applied across 
the longitudinal and circumferential strain receiver 
windings connected in series. The electrical output of 
the strain receiver system, taken from the common lead 
of the receiver, is applied as the input to the coupling 
system amplifier. The output from this amplifier is the 
coupling system output that acts as the input signal to 
the indicating system. 
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OSCILLOSCOPE 








PHOTOGRAPH SHOWING PARTS OF COUPLING SYSTEM EXTERWAL TO. 
RECEIVER UNIT AND CATHODE RAY OSCILLOSCOPE AS INDICATING SYSTEM 


Fig. 18 is a photograph showing a cathode-ray oscil- 
loscope, which is used as the indicating system, and the 
box containing coupling system parts that are outside 
the pressure receiver unit. 


Mechanical Features of Pressure Receiver Unit 


The upper part of the photograph of Fig. 19 shows the | 
external appearance of a complete catenary diaphragm- 
strain generating tube pressure receiver unit; the lower | 
part shows the four component parts of the pressure re- 
ceiver unit. Figs. 20 and 21 give cross-sectional dia- 
grams showing the operating relationships among in- 
ternal parts of the pressure receiver unit. Correspond- 
ing parts are given the same labels in Figs. 19, 20, and 
21. The catenary diaphragm is formed as a separate 
unit and is attached to the shell by means of a forcing 
fit on a silver solder joint. The shell carries an outside 
thread to fit into an 18-mm. spark plug hole, a locating 
shoulder to locate the cooling air guide tube, and an in- 
ternal thread to receive the connector sleeve that 
serves to hold the component parts of the unit together. 
The cooling-air guide tube has a flange that fits against 
the locating shoulder. When the assembly is complete, f 
the cooling-air guide tube forces air, supplied under 
pressure, through the connector sleeve to flow partly 
along the length of the strain generating tube and partly 
outward across the catenary diaphragm through cool: | 
ing-air passage holes in the strain generating tube be- 
fore leaving to the atmosphere through the cooling-aif 
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FIG 20. CROSS SECTION OF CATENARY DIAPHRAGM PRESSURE RECEIVER UNIT 


exit holes in the shell. The path of this air circulation is 
shown in Fig. 21. The strain generating tube support 
fits into the shell on top of the cooling-air guide tube 
flange. A tube position adjusting spacer is used to lo- 
cate the strain generating tube support flange so that 
the rounded end of the strain generating tube, which acts 
as the force summing member, is brought into proper 
contact with the catenary diaphragm. In practice, 
the tube is located so that the diaphragm is deflected 
from its equilibrium position by approximately 0.002 
in. The strain receiver winding is cemented to the out- 
side of the strain generating tube with the arrangement 
shown in Fig. 15. Wire connections attach the ends 
and the midpoint of the strain receiver winding to a 
three-prong connector, which is held in place by a shoul- 
der in the connector sleeve when this sleeve is screwed 
into the shell. The complete pressure receiver unit is 
approximately 2 in. long and will fit into holes designed 
for 18-mm. spark plugs. 
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(VI) PERFORMANCE REALIZED BY AN INDICATOR WITH 
CATENARY DIAPHRAGM PRESSURE RECEIVER UNIT 


Sensitivity of the Pressure Receiver Unit 


An indicator of the type described in Figs. 19 and 20 
can be operated to receive pressure as its actuating in- 
put and constant voltage as its modifying input and to 
Calibration 
data taken on a number of catenary diaphragm indica- 
tors designed to operate under 1,000 lbs. per sq.in. have 
shown that the realized modified voltage input, pressure 
input - voltage output sensitivity is about 0.001 millivolt 
per Ib. per sq.in. Over the designed operating range of 
10 lbs. per sq.in. below atmospheric pressure to 1,000 
Ibs. per sq.in. above atmospheric pressure, the devia- 


produce voltage as its essential output. 


tion of this sensitivity from its specified value is always 
less than 1 per cent. 


Comparison of Catenary Diaphragm Unit with Flat- 


Diaphragm Unit 

The operating characteristics achieved by the use of 
the pressure receiver design features described in Figs. 
19 and 20, as compared with the results from a flat- 
diaphragm pressure receiver, were measured by taking 
records from units of the two types under identical 
conditions. The flat-diaphragm unit used for compari- 
son purposes was a capacitance-type indicator developed 
by General Motors Research Corporation.? This unit, 
which represents the best available equipment of its 
type, is referred to as the flat-diaphragm pressure re- 
ceiver in the figures that follow. 
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FIG. 21 CROSS SECTION OF CATENARY DIAPHRAGM PRESSURE RECEIVER UNIT 
SHOWING PATH OF COOLING AIR. 
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VIBRATION RESPONSE PRE RE INDICATOR WITH CATENARY O/4PHR AGM 

PR RE RE 8 
FIG 22 COMPARISON OF VIBRATION RESPONSE OF CATENARY DIAPHRAGM 
PRESSURE RECEIVER WITH VIBRATION RESPONSE OF FLAT DIAPHRAGM 


PRESSURE RECEIVER 


Acceleration sensitivity effects were tested by subject- 
ing the support carrying the comparison pressure re- 
ceiver units to shocks produced by hammer blows. Os- 
cillograph records from the comparison units, obtained 
with shocks of this kind, are shown in Fig. 22. The 
upper record in the figure is from the indicator using the 
flat-diaphragm indicator; the lower record is from the 
catenary diaphragm unit indicator. These records 
show that the response of the catenary diaphragm unit 
is about one-fifth as much as the response of the flat dia- 
phragm unit. This result is in accord with theoretical 
predictions based on the sensitivity ratio curves of Fig. 
13. 

The effect of temperature gradient through the dia- 
phragm material was demonstrated by subjecting the 
pressure receivers of the comparison indicators to the 
same fluctuation in temperature by repeatedly covering 
and uncovering a gas flame playing on the pressure 
seals. The curves of Fig. 23 show that the effect of a 
rapidly changing temperature on the catenary dia- 
phragm is less than one-twentieth as much as the effect 
of the same temperature change on the flat diaphragm. 
This performance improvement is substantial but is less 
than that predicted from theory for a perfectly flexible 
pressure seal diaphragm. The difference is due to the 
deviation of the actual catenary diaphragm and strain 
generating tube system from an ideal system. 

The effect of average temperature gradient in chang- 
ing the zero reference level in terms of indicated pres- 
sure is shown in the curves of Fig. 24. It appears from 
these curves that, for the indicator with the flat dia- 
phragm pressure receiver (with water cooling), the zero 
reference point for the pressure indications shifted by an 
amount almost equal to the full-scale reading when the 
temperature of the support was changed approximately 
200°F. On the other hand, under the same temperature 
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change, the indicator with the catenary diaphragm 
pressure receiver showed a pressure reference level 
change of approximately 1 per cent of that for the flat- 
diaphragm unit. This result is in accord with the theo- 
retical prediction that, with a diaphragm of high flexi- 
bility, average temperature changes would not affect the 
indicated pressure output. 

The test results of the effect of average support tem- 
perature on the pressure input - voltage output sensitivity 
of the catenary diaphragm pressure receiver unit is 
given in Fig. 16, together with the theoretical calcula- 
tions. The pressure - voltage sensitivity performance 
ratio, which is the ratio of the actual sensitivity to the 
sensitivity for the unit at a temperatureof 68°F.,showsa 
deviation of less than 1 per cent from unity for tempera- 
tures between 68° and 400°F. With air cooling, the 
operating temperature of a catenary diaphragm pres- 
sure receiver never exceeds about 150°F.; therefore, 
the system sensitivities of indicators using catenary 
diaphragm pressure receiving units will be free from 





FIG.23. RESPONSE F PRE RE INDICATORS TO PERIODICALLY VARYING TEMPERATURE 


SHIFT, IN psi 














SUPPORT TEMPERATURE, IN °F 


TEMPERATURE ON ZERO REFERENCE OF PRESSURE 
MAXIMUM SCALE READING 


FIG. 24 EFFECT OF SUPPORT 


NDICATORS WITH 1000 ps 








ee eae 





Stan Re 


hragm 
> level 
1e flat- 
> theo- 
1 flexi- 
ect the 


t tem- 
sitivity 
init is 
ilcula- 
mance 
to the 
10WSa 
pera- 
g, the 

pres- 
-efore, 
enary 

from 


| 





‘TURE 


~ 300 



















A NEW 





HIGH-PERFORMANCE ENGINE 











607 





INDICATOR 








CATENARY 
DIAPHRAGM rac BRIDGE ce, satel RC aan ne 
' ry 


——  0Ul ll Kl OS OO 





DEMODULATING ciRcUIT FILTER 


CARRIER 


— Se REJECTION] (-D-C AMPLIFIER] FOSCILLOSCOPE] 


4 








GAGE UNIT 











_—_—+ — 











SUBSEQUENT STAGES OF 
R-C AMPLIFIER 



































OSCILLATOR 
L FOR &C —) 
EXCITATION 


FIG. 25 


temperature effects under all ordinary operating condi- 
tions. 

Life tests of 100-hours engine operation have shown 
no essential deterioration of the units. 


(VII) TypicaAL RECORDS OF ENGINE PRESSURE 
VARIATIONS TAKEN BY INDICATORS USING THE 
CATENARY DIAPHRAGM PRESSURE RECEIVER 


Pressure - Time Records 

Pressure records from operating engines may be 
taken by applying to the receiver winding of a catenary 
diaphragm pressure receiver either a constant voltage, 
as shown in the wiring diagram of Fig. 17, or an a.c. 
voltage, as shown in the wiring diagram of Fig. 25. A 
typical pressure - time record is produced when constant 
voltage excitation is used and when the sweep frequency 
of a conventional cathode-ray oscilloscope, acting as the 
time axis, is properly synchronized with the cycle fre- 
quency. A record of this kind is visible on the face of 
the cathode-ray tube pictured in Fig. 18. In taking 
permanent photographic records, it is more accurate to 
provide the time axis by means of a constant speed film 
rather than by using the sweeping circuit of the oscillo- 
scope. The overall accuracy of single-cycle pressure 
records taken by photographing the cathode-ray trace 
can be checked by comparing points taken from these 
records with corresponding points on an average pres- 
sure record taken at the same time as the continuous 
pressure records. Fig. 26 shows points from a typical 
single-cycle pressure record, taken by an indicator us- 
ing the catenary diaphragm pressure receiver with con- 
stant voltage excitation, superimposed on the corre- 
sponding average pressure record from an M.I.T. Bal- 
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THE PHASE SENSITIVE DEMODULATING CIRCUIT OF THIS 
FIGURE WAS ADAPTED FROM AN UNPUBLISHED WORK 
OF PROF. R.MUELLER OF THE INSTRUMENTATION 
LABORATORY, M.I.T. 


CIRCUIT DIAGRAM OF A-C CARRIER SIGNAL GENERATING SYSTEM FOR CATENARY DIAPHRAGM PRESSURE INDICATOR. 


anced Pressure Diaphragm Indicator.** This figure 
shows that all the points from the single-cycle record 
fall within the spread of averaging indicator points pro- 
duced by the inherent variations between successive 
engine cycles. 

This agreement has been consistently repeated for 
variations of engine conditions. In general, it appears 
that the continuous pressure records for an indicator 
using the catenary diaphragm pressure receiver will 
have inaccuracy components no greater than the inac- 
curacies of high-performance averaging indicators. 

Constant voltage excitation for the strain receiver 
winding is easily obtained in practice by using conven- 
tional ‘“‘B”’ batteries. If the frequency of the engine 
power cycles were faster than 4 cycles per sec., the char- 
acteristics of amplifiers usually included in commercial 
cathode-ray oscilloscopes would give satisfactory pres- 
sure records. However, at engine speeds corresponding 
to less than 4 power cycles per sec., a satisfactory sys- 
tem requires amplifier performance not usually built 
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FIG 27 ENGINE PRESSURE - TIME RECORD FROM INDICATOR USING CATENARY DIAPHRAGY 
PRESSURE RECEIVER UNIT WITH ALTERNATING GURRENT EXCITATION AT 5000 CYCLES 
PER SEGOND. 
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into standard equipment. The practical result of this 
defect in indicating system performance is that the pat- 
tern on the cathode-ray tube is considerably distorted, 
with insufficient screen indication during lower rate of 
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pressure change. This effect may be eliminated by us. 
ing a relatively high-frequency (5,000 cycles per sec. or 
higher) a.c. voltage for exciting the strain receiver wind- 
ing. When a carrier frequency of this kind is used, the 
pressure modulated output of the pressure receiver unit 
can be transmitted by commercial resistance-capaci- 
tance coupled amplifiers, without causing the indicated 
pressure pattern to shift its zero when changes occur in 


engine-operating conditions. By incorporating a de- 


modulating circuit in the electrical components of the : 


coupling system or the indicating system, as shown in 
Fig. 25, the pressure indications will appear as a simple 
trace on the face of the cathode-ray tube. It is also pos- 
sible to apply the modulated signal output from the 
strain receiver winding directly to the input terminals of 
a cathode-ray oscilloscope indicating system. Because 
of the lower linear velocity of the cathode-ray spot 
near the peak of its travel, the envelope of the pattern 
is outlined by bright points, corresponding to longer 
times for the electron beam to excite the fluorescent 
screen. Fig. 27 is a photograph showing the pressure - 
time record produced by alternating-current excitation 
of the strain receiver winding without use of demodula- 
tion. Records of this kind can be generated by a com- 
paratively simple electrical coupling system, with a com- 
mercial cathode-ray oscilloscope for the indicating sys- 


tem. 


Pressure - Volume Records 


Pressure - volume records can be obtained from a pres- 
sure indicator system like that of Fig. 18 by making the 
horizontal motion of the cathode-ray spot follow the 
piston motion instead of being a linear function of time. 
A simple converter from crank angle into signal of pis- 
ton position can be made by using a model crank and con- 
necting rod combination, with relative proportions iden- 
tical with those of the actual engine, to stretch a spiral 
spring in proportion to piston displacement. If the 
crank of the model is made to rotate with the crank of 
the engine, the electrical signal output from the strain 
of the stretching spring will be proportional to the en- 
gine piston position. A pressure-volume record like 
that of Fig. 28 is produced when this piston position 
signal is applied to the horizontal deflection plates of 
the cathode-ray oscilloscope, whose vertical deflection 
plates receive the signals from the pressure receiver. 


Detonation Records 

The sample records of Figs. 26, 27, and 28 illustrate 
the ability of indicators using the catenary diaphragm 
pressure receiver to measure the pressure variations as- 
sociated with the power-generating functions of an en- 
gine. Any generally satisfactory continuous pressure 
indicator must not only record average cycle variations 
but must also accurately measure the relatively high-fre- 
quency fluctuations associated with detonation. Fig. 
29 shows a pressure - time record taken from an engine 
operating with severe detonation. The time scale is ex- 
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panded longer than that for power analysis in order to 
show more clearly the fine structure of the high-fre- 
quency pressure changes. Amplitudes and frequen- 
cies of the pressure variations associated with detona- 
tion are marked in Fig. 29 as functions of crank 
angle. 

A curve showing detonation frequency as a function 
of crank angle, obtained from the record of Fig. 29, is 
given in Fig. 30. Also plotted in’ Fig. 30 are curves of 
detonation frequency as a function of crank angle, based 
on theoretical computations* and on data taken from 
records made with the same engine type and operating 
conditions from an indicator using a flat-diaphragm 
pressure receiver. * 

The detonation frequency curve, based on data from 
the indicator using the catenary diaphragm, matches 
the computed results and also matches the measurement 
made by an indicator differing greatly in design from 
the catenary diaphragm unit. In addition to the agree- 
ment in frequency, the detonation pressure amplitudes 
of Fig. 29 have magnitudes that are approximately 
equal to estimated amplitudes derived from measure- 
ments of rates of change of pressure.*® For example, 
these measurements showed that, under detonating 
conditions similar to those used in making the record 
of Fig. 29, the detonation pressure amplitude at 20° af- 
ter top dead center was about 40 Ibs. per sq.in. Refer- 
ence to Fig. 29 shows that this amplitude is substan- 
tially equal to the directly measured amplitude for the 
same crank angle. 


SUMMARY 


It has been brought out in this paper that continuous 
pressure indicators available up to this time have per- 
formance defects that make them unsuitable for obtain- 
ing satisfactory pressure records from modern engines. 
These difficulties have been largely due to undesirable 
characteristics inherent in the flat-diaphragm type of 
The discus- 
sion of this paper summarizes the problem to show that 
excessive sensitivity to temperature, mechanical strain 
resulting from mounting stresses, and mechanical vi- 
brations in the walls of the pressure chamber are the 
principal sources of trouble. A semiquantitative discus- 
sion leads to the conclusion that the flat-diaphragm 
type pressure receiver will always be subject to intoler- 
able inaccuracies under the operating conditions of 


pressure receiver used with the indicators. 


modern engines. 

A new pressure receiver with greatly improved per- 
formance is described. This instrument is practically 
free from temperature and installation-strain effects 
and has a greatly reduced sensitivity to vibration as 
compared with the flat-diaphragm type pressure receiv- 
ers. This unit is made up of a catenary diaphragm and 
a strain generating tube, wound with a wire strain gage. 
Typical results from laboratory tests and engine runs 
are included to show that, on the basis of an inaccuracy 
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tolerance of 1 per cent, indicators using the new receiver 
will be satisfactory under engine operating condi- 
tions. 
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Some Aspects of Convertible Aircraft Design 


ROBERT L. LICHTEN* 
Bell Aircraft Corporation 


SUMMARY 


This paper is a contribution toward a quantitative engineering 
treatment of convertible aircraft design. A brief discussion of 
design requirements and proposed configurations is given. The 
inalysis covers the following: 

1) Conversion from Helicopter Rotor to Propeller Operation.— 
A semigraphical process for computing rotor characteristics 
during the conversion process when blade-section inflow angles 
Calculations 
for a typical case indicate that no unusual or undesirable rotor 
behavior occurs, blade flapping decreasing progressively from 
the magnitude existing in helicopter operation. The rotor is 
shown to operate efficiently throughout, indicating that aircraft 
performance is not impaired during the conversion interval. 

2) Rotor Aerodynamic Characteristics in Propeller Operation.— 
It is shown that satisfactory propulsive efficiency can be ob- 
tained with isolated rotors having relatively low solidity and 
blade twist compared to conventional propellers. In order to 
obtain good efficiency at high-forward speeds, it is necessary to 
operate the rotor at high pitch and power coefficient. In prac- 
tice, this means that when a helicopter rotor is converted to 
propeller operation its tip speed must be reduced considerably. 
Model propeller test data are compared with calculations based 
on modified vortex theory, and it is shown that the calculation 


can no longer be considered small is presented. 


method yields reasonably accurate results for the somewhat un- 
usual operating conditions of the rotor-propeller. 

(3) Effect of Gravity on a Lag-Hinged Rotor in Propeller Opera- 
tion.—The undamped in-plane motion due to gravity excitation 
of a lag-hinged rotor blade is analyzed for operation with rotor 
axis horizontal. In the case of a three-bladed rotor, it is shown 
In the 
case of a two-bladed rotor, the amplitude of the vibratory force 
applied to the rotor hub is considerably greater than total blade 
weight, and therefore a two-bladed convertible rotor probably 
must be of the semirigid or rigid type to avoid this condition 


that the resulting hub vibratory loads are negligible. 


INTRODUCTION 


ne AIRCRAFT, which combine the ability of 
the helicopter to take off vertically and hover with 
the ability of the airplane to fly at high speed over long 
distances, have been proposed in various forms for 
many Design 
high-speed flight have been generally known for some 
time; 


years. requirements for moderately 
however, it was not until the more pressing 
problems associated with the development of practical 
helicopters were solved and corresponding design knowl- 
edge for the low-speed part of the flight spectrum be- 
came available that such proposals could be given seri- 
ous consideration. 

In the writer's opinion, the primary requirements for 
a successful convertible aircraft are that low- and high- 
speed performance be essentially as good as the per- 
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formance of comparable helicopters and airplanes in 
their respective flight régimes and that the excellent 
controllability and maneuverability of present-day 
helicopters at zero and low air speeds be retained. 
As an airplane, the convertible aircraft must be able 
to fulfill an economic function. At low speed, it 
must be able to execute rapid and precise changes in 
air speed and air position, must be able to make safe 
landings after power failure in ground areas as small 
and restricted as those normally used for power-on 
operation, and must be relatively undisturbed by gusty 
air. Failure to meet these requirements would result in 
a degree of safety and utility too low to justify the cost 
of development. 


DISCUSSION OF TYPES 


A brief discussion of the more prominent proposed 
types of convertible aircraft is given here to indicate 
the reasons for the direction taken in the analysis. It 
should be emphasized that all statements made are 
matters of opinion at this time. 

Two basic classes of direct-lift aircraft may be de- 
fined: (A) aircraft in which direct lift is obtained by 
imparting a relatively large downward acceleration to 
a relatively small mass flow of gases, as by use of 
rocket or turbojet engines; (B) aircraft in which direct 
lift is obtained by imparting a relatively small down- 
ward acceleration to a relatively large mass flow of at- 
mospheric air,’as by use of helicopter or cyclogyro rotor 
systems. 

Class (A) fails to meet the previously stated require- 
ments for convertible aircraft in that power failure 
during direct-lift operation would be disastrous, and 
therefore is not considered further. Class (B) will be 
limited to use of helicopter-type rotor systems for this 
discussion, and convertible aircraft in this class may be 
subclassified into two main groups, each having two 
subgroups, as follows: 

(1) Aircraft having one or more helicopter-type 
rotors as the principal source of lift in low-speed flight; 
for high-speed flight, rotor rotation is stopped. (a) 
Longitudinal fuselage axis is approximately vertical at 
take-off and approximately horizontal at high speed; 
(b) longitudinal fuselage axis is approximately hori- 
zontal throughout the entire flight range. 

(2) Aircraft having one or more _ helicopter-type 
rotors with axes approximately vertical as the prin- 
cipal source of lift in low-speed flight; for high-speed 
flight, rotor axes are tilted to an approximately hori- 
zontal attitude in which the rotors act as propellers. 
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(a) Longitudinal fuselage axis is approximately vertical 
at take-off and approximately horizontal at high speed, 
and rotor axes are fixed with respect to fuselage; (b) 
Longitudinal fuselage axis is approximately horizontal 
throughout the entire flight range, and rotor axes tilt 
with respect to fuselage. 

Types (la) and (2a) represent entirely possible 
solutions. Type (la) has been embodied in proposed 
aircraft having ram-jet or pulse-jet engines mounted 
at the tips of relatively rigid blades to drive the rotor 
for helicopter operation; for high-speed flight the 
blades are feathered and act as wings, while the jet 
engines provide direct propulsion. Type (2a) has 
been embodied in many designs, such as those pro- 
posed by Zimmerman, Young, Leonard, and Robins. 
In both these types, conversion from helicopter to air- 
plane flight and return may be expected to require either 
a discontinuous maneuver, such as a dive or zoom, 
or a continuous maneuver in which control and excess 
power available are likely to be marginal for an ap- 
preciable time interval. The required 90° range of 
fuselage attitude would probably make such an air- 
craft unsuitable for carrying more than one or two 
persons. Autorotational descent in the event of power 
failure, with sufficient forward speed to minimize 
sinking speed and to provide a choice of landing sites, 
would be a difficult and awkward maneuver at best. 
For type (2a), the coaxial rotor arrangement, with two 
superimposed oppositely rotating rotors, is most often 
used, leading to solidity values that are objectionably 
high for propeller operation, as will be shown later. 

An early example of type (1b) was the Herrick Con- 
vertaplane demonstrated with partial success in 1937. 
Other designs of this type have been proposed by Wil- 
ford and Gazda. Such aircraft generally require sepa- 
rate power transmission systems for driving the rotor in 
helicopter flight and for propelling the aircraft after 
stopping the rotor. Simplification in this respect might 
be effected by use of some form of jet drive for both 
configurations. A one-bladed rotor offers advan- 
tages for this type; however, such a rotor presents 
many new problems even for use on conventional heli- 
copters. Mechanical requirements for stopping and 
starting large rotors in flight are likely to be severe. 

Type (2b) embodies, in general, the most conven- 
tional configurations as suggested by Blount, LePage, 
Stuart, DuMonge, and others. This type appears to be 
the most suitable for large cargo and personnel trans- 
port aircraft but has the disadvantage of a relatively 
heavy structure and power transmission system. How- 
ever, it introduces the least new design problems, per- 
mits adequate provision for safe operation after power 
failure, and seems most promising for immediate de- 
velopment. Therefore the analysis presented in the 
following sections is directed primarily toward this 
type, although it is applicable in many respects to other 
types discussed above. 

Additional discussion is given in reference 9. 


The plane of the rotor disc is a plane perpendicular to the rotor 
shaft axis passing through the rotor hub. 
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= tan! (Cdg/ C1) 


‘rotor shaft power coefficient = P/[(1/2)e(wR)*xR?] 
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NOTATION 


slope of blade section lift coefficient curve, per rad 

coefficient of cos my term in expression for 8 

coefficient of cos my term in expression for 4 

number of blades per rotor 

coefficient of sin my term in expression for 3 

tip-loss factor 

coefficient of sin my term in expression for y 

chord of constant-chord blade, ft. 

blade section profile drag coefficient 

blade section lift coefficient 

rotor thrust coefficient = 7/![(1/2)p(wR)?xR?] 

rotor diameter, ft. 

distance from rotor shaft axis to lag hinge, ft. 

Goldstein propeller theory factor («x in reference 2 

centrifugal force on one blade, lbs. 

force on rotor hub in positive y direction, Ibs. 

force on rotor hub in positive z direction, lbs. 

moment of inertia of blade about lag hinge, slug-ft.? 
W,R? ath 

y= =e 

g é: + Le 





WLR /2w? 

distance from lag hinge to blade center of gravity, | 

ft. : 

distance from lag hinge to point of application of § 
resultant Fy, ft. 

rotor rotational speed, revolutions per sec. 

rotor shaft power, ft.lbs. per sec. 

pitching angular velocity of rotor shaft axis, positive 
nose up, rad. per sec. 

radius to a blade element, ft. 

radius to blade tip, ft. 
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= rotor thrust, vector coincident with rotor shaft axis, ' 


= true air speed along flight path, ft. per sec. 
= total 


lbs. 

component of resultant velocity at a blade element, 
parallel to rotor shaft axis, ft. per sec. 

component of resultant velocity at a blade element, 
normal to blade span axis and parallel to rotor disc 
plane, ft. per sec. 

average induced velocity 
rotor disc plane, ft. per sec. 
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at rotor disc, normal to 


induced velocity at a blade element (profile 
drag neglected), ft. per sec. 

resultant of V and wr at a blade element, ft. per sec. 

resultant of IV) and w at a blade element, ft. per sec. 

weight of one blade, Ibs. 

r/R 

fixed horizontal coordinate in plane of rotor dise, 
propeller operation, ft. 

fixed vertical coordinate in plane of rotor disc, propeller 


Se EE 


operation, ft. 

angle between flight path and rotor disc plane, posi 
tive when leading edge of disc is raised, rad. 

blade element angle of attack relative to zero-lift 
chord line, rad. 

flapping angle between blade span axis and rotor dist 
plane, positive upward, rad. ; 

lag angle between projections in rotor disc plane o 
blade span axis and line connecting rotor shaft axis 
and lag-hinge, rad. 

tan! (w/W), rad. 

propulsive efficiency, propeller operation, = TV/P 

rotor profile efficiency = T[v — Vsin (a + a;)]/P 
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r) = pitch angle between blade section zero-lift chord 
line and rotor disc plane = 0 + 6x, rad. 
d = inflow ratio = (V sina — v)/wR 


ll 


advance ratio = (V cos a)/wR 
air density, slugs per cu.ft. 


p = 

e = rotor solidity = bc/rR 

do = tan~! (V/wR), rad. 

@ = go + ¢, rad. 

y = blade azimuth angle measured from downwind or top 
position, rad. 

w = rotor rotational speed = dy /dt, rad. per sec. 


LIFTING ROTOR TO PROPELLER 
OPERATION 


CONVERSION FROM 


Introduction 


In converting from lifting rotor to propeller opera- 
tion in flight, it is highly desirable that the transition 
be made without requiring unusual piloting technique, 
penalizing performance appreciably, or causing erratic 
The type of maneuver envisioned for 
the conversion process is as follows: The aircraft at- 
tains maximum forward speed as a helicopter. The 
pilot then actuates the conversion mechanism, causing 
the rotor axis to swing from an approximately vertical 
to an approximately horizontal attitude parallel to 
the flight path. This process occupies a time interval 
of the order of 15 sec., at the end of which the rotor 
is operating as a propeller and the aircraft weight load 
is carried by fixed wings. Conversion from propeller 
to lifting rotor operation is the reverse of the process 
described above. 


rotor behavior. 


Experimental Data 


Flight tests of a convertible aircraft model are re- 
ported in reference 7. Using the classifications given 
previously, this model was of type (2a). It had a 39-in. 
diameter, two-bladed rotor, and gimbal mounted to 
permit flapping and feathering and was driven by an 
electric drill motor. A 48-in. span by 8-in. chord 
wing was mounted below the rotor with chord line 
parallel to rotor shaft axis. The model was mounted 
at the end of a dural ‘‘fish pole’ 18 ft. long and was 
flown in a circular path around the operator. The 
model was free to pivot about the pole axis, and the 
operator had longitudinal cyclic feathering control, 
collective pitch control, and motor torque control. 
Thus, the operator controlled the model about the 
pitch axis by means of rotor thrust vector inclination, 
but the model was restrained from yawing and ‘rolling. 

In flight, the transition from hovering to high-speed 
flight with rotor axis horizontal could be made rapidly 
or slowly and was smooth and continuous. The model 
could be flown steadily with the rotor shaft axis in a 
range of attitudes from vertical to horizontal. It was 
also possible to perform loops and vertical figure eights. 
Rotor behavior was normal throughout all maneuvers. 

Approximate measurements of power input gave 
evidence of a sharp increase in power required in the 
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from model flight-test data. Model wing chord line fixed paral- 


lel to rotor axis. 


vicinity of a = —60° over that required in hovering 
and in forward flight at smaller values of a. This is 
attributed to operation of the wing in a stalled attitude 
at the relatively high forward speed corresponding to 
this axis inclination. As axis inclination is further 
increased, the wing would be expected to unstall with a 
consequent drop in power required, followed by a rise 
as maximum forward speed is attained with a> —90°. 
Such a trend is exhibited by the test data. A curve 
showing the approximate variation of power required 
with rotor axis inclination at one rotational speed is 
This curve was cross-plotted from 
No air- 


given in Fig. 1. 
faired curves based on the reference 7 data. 
speed data were recorded. 


Analysis 


During the conversion process, the principal condi- 
tions of operation which might be expected to lead to 
unusual rotor behavior because of differences from 
steady-state lifting rotor or propeller operating condi- 
tions are: (a) airflow interference caused by nonrotor 
parts of the aircraft; (b) the pitching angular motion 
of the rotor hub as the rotor axis is converted from an 
approximately vertical to an approximately horizontal 
attitude; (c) the large inflow angles, together with con- 
tinuous variation of inflow angle around the blade path, 
which are experienced by the blade sections. 

Interference effects depend on the particular con- 
figuration of the subject aircraft. They may be con- 
sidered as superimposing on the airflow pattern, sur- 
rounding an isolated rotor in a given operating condi- 
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tion, a velocity and static pressure variation at each 
point in the pattern caused by the influence of other 
parts of the aircraft. Good design will minimize such 
effects; however, there are instances of difficulties en- 
countered with propeller installations in conventional 
aircraft due to interference. In the present case, rotor 
behavior difficulties due to interference are not ex- 
pected to appear, since such a rotor is already designed 
to accommodate the large variations with blade azi- 
muth position and forward speed of section angle of 
attack and relative airflow velocity experienced in 
helicopter operation. Therefore interferencé effects 
are not herein considered further, although they may 
have a considerable effect on stability and control in 
particular cases. 

It may be shown, exactly in the case of a flapping- 
or teetering-type lifting rotor at zero air speed and 
approximately for such a rotor in forward flight, that a 
steady pitching angular velocity of the hub will pro- 
duce lateral flapping 


—P 


b = — pw/w = 


due to aerodynamic damping of the pitching motion 
and will produce longitudinal flapping due to blade 
mass whose magnitude is dependent on the ratio of 
blade moment of inertia about the flapping hinge to air 
density and which, for typical rotor blades, is approxi- 


mately 


These values should still be good approximations for 
flapping due to hub pitching velocity during the con- 
version process. For an angular motion of 90° in 15 
sec., about 0.10 rad. per sec., and a typical rotor speed 
w = 25 rad. per sec., the resultant flapping amplitude 


would be 


Va? + b,2 = 2.7p = 0.011 rad. = 0.6° 


This may be considered negligible so far as any sig- 
nificant effect on rotor behavior is concerned, and this 
conclusion is supported by the fact that present small 
helicopters can obtain rates of roll or pitch in excess of 
| rad. per sec. in maneuvers without experiencing ab- 
normal rotor behavior. Therefore it has been con- 
sidered permissible to neglect the effect of hub pitch- 
ing velocity in analyzing rotor behavior during con- 
version, which permits a quasi-static treatment in 
that steady-state operation at rotor angles of attack 
intermediate between lifting rotor and propeller condi- 
tions may be assumed. 

The aerodynamic analysis of rotor operation in this 
intermediate range is complicated because the inflow 
angle cannot be considered a small angle as in lifting 
rotor analysis, while conditions at any blade station 
do not remain constant with respect to blade azimuth as 
in propeller analysis. This leads to a strip analysis 
with graphical integration for rotor characteristics. 
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Schematic diagram of a typical blade section for con- 
version analysis. 


Fic. 2. 


In addition to the assumptions that interference ef- 
fects and hub pitching velocity may be neglected, the 
analysis is based on the following: (a) Only the com- 
ponent of induced velocity normal to the rotor disc 
plane need be considered. For this component a mean 
value over the disc as given by momentum theory may 
be used. (b) only the constant and first harmonic 
components of the blade flapping motion need be con- 


sidered. (c) Eq. (2) of reference 8, 


tan a = d/p + Cr/4B*p (A? + yw?) (I 


remains valid during conversion. This equation in- 
cludes an approximation for lifting rotor induced flow, 
based on simple momentum theory and analogy with 
finite airfoil theory. 

The above assumptions are equivalent to those used 
in standard blade-element analysis of lifting rotors. 
Accordingly, 

Wm =x+usiny (2) 
Up = X\ — x(q sin Y — Bb cosy) — 
ucos y (do — a, cosy — h sin y) (3) 


Referring to Fig. 2, 


a = A “f- Ax + tan 


1 (uy/U;) 


l dCr _ - ° 
= (u?? + uy”) (¢, cos 6 + Ca sin @) (4) 
o dx 
Pe Oe Sani , 4 
— = x(u;? + Up”) (Ca cos d — C, Sin d) (5) 
o dx 


In applying these expressions, blade airfoil section 
characteristics are taken from a plot that covers a large 
range of a, extending well into the negative stall re- 
gion. 

No further details of the analysis will be given ex- 
cept for the following outline of the procedure followed: 

(a) For a given rotor, choose a forward speed and 
value of a for analysis. A value of u is derived. 

(b) Select a range of values of \ and % appropriate 
to the chosen a. 

(c) Assuming zero blade flapping, calculate Cy, for 
each possible combination of A and 4 by evaluating 


Eq. (4) over a range of values of x at intervals of ¥ 
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around the blade path (every 30° is recommended). 
The resulting values of dC;/dx at each value of x 
are averaged with respect to y, and C, is obtained 
by graphical integration. The tentative assumption 
is made that thus-obtained values of Cy, are lit- 
tle affected by introduction of actual blade flap- 
ping. 

(d) Sets of A, wu, and C; values from (c) are substi- 
tuted in Eq. (1), yielding values of a which are plotted 
versus the assumed A and @ values and the calculated 
Cr values. Intersections at the originally assumed 
a value determine possible operating conditions (sat- 
istying both momentum and blade-element condi- 
tions). 

(e) For each of the possible conditions, a substitute 
effective value of \ is determined by substituting the 
corresponding Cp and @ in Eq. (14) of reference 8. 
This is then used to compute flapping coefficients do, 
a,, and 6, from Eqs. (11) to (13) of reference 8. 

(f) Using these flapping coefficients, values of C; and 
Cp are computed following the procedure outlined in 
(c) for the range of possible operating conditions. 
Using the dC,/dx values obtained during this compu- 
tation, values of x(dC;/dx) are computed and integrated 
over the blade radius for each blade azimuth position. 
This quantity is proportional to the aerodynamic thrust 
moment about the flapping or teetering hinge, and a 
harmonic analysis of its variation with y is performed. 
A corresponding procedure is followed for the zero- 
flapping cases. The correct values of a; and }; are 
determined by interpolation between the two sets of 
first harmonics obtained from the 


Sc? x(dCy/dx)dx 


variation with y. It may be assumed that the sine 
component varies linearly with a; and the cosine com- 
ponent linearly with },; for a hinged rotor, a; and }, 
are such that the first harmonic components are re- 
duced to zero. In the actual calculations made, it was 
found that the flapping coefficients obtained in step (e) 
produced this result directly within the limits of ac- 
curacy of the method so that interpolation was unneces- 
sary. This is indicated by the example shown in Fig. 3. 
In all cases the final flapping values should be suf- 
ficiently close to those calculated in step (e) that the Cr 
values obtained in step (f) do not require recalculation. 
This conclusion is supported by the small changes be- 
tween values for zero and for calculated flapping shown 
in Table 1 for the example case. It has been assumed 
that similar reasoning may be applied to Cp. 


TABLE 1 





a = —60°, » = 0.15, 0 = —20° 
4 (rad.) 0.600 0.630 0.660 
A —0.2600 —0.2613 — 0.2627 
Cr/e (zero flapping) —0.0125 0.0328 0.0822 


Cr/o (calculated flapping) | —0.0161_ 0.0334 0.0827 
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Fic. 3. Calculated variation of thrust moment coefficient with 
azimuth angle for example rotor. a = —60°; uw = 0.15; & = 
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coefficients with blade pitch for example rotor. a = —60°; 
nw = 0.15; o = 0.04; 6 = —20°. 
(g) In the a = —90° case (propeller operation), a 


considerable simplification can be made, since there is no 
variation with y and first harmonic flapping is zero. 
The general method outlined above may still be fol- 
lowed in order to obtain a comparative set of results 
for this limiting condition. 

Some results for a set of calculations with a = —60° 
are presented as anexample. Table 1 shows a compari- 
son of C, values as computed for cases with zero and 
with calculated flapping. Differences are within the 
limits of accuracy of the computations. 

In Fig. 3, the variation of thrust moment coefficient 
about the flapping hinge is shown for these cases. It 
may be noted by inspection that the calculated flap- 
ping reduces the first harmonic variations essentially to 
zero. In Fig. 4, final rotor characteristics for this 
operating condition are shown. 
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TABLE 2 
Rotor disc loading = 3.5 lbs. per sq.ft. (in helicopter operation) 


Rotor power loading = 


9.8 lbs. per shaft hp. 
.04 


wR = 612 ft. persec., o = 








p = 0.00238 slug per cu.ft., 0, = —20 

a V 0 a i 7 

(Deg.) ‘Flight Condition (M.P.H.) (Deg.) ” nN (Deg.) (%) ” 
0.0 Hovering 0 27 .0 0.0 —0.0665 0.0 100.0 0.72 
— 20.3 Helicopter high speed 133 28.8 0.30 —0.1130 8.05 100.0 0.77 
—30.0 Conversion 125 30.0 0.26 —0.1542 6.51 61.4 0.83 
—60.0 Conversion 125 36.8 0.15 —0.2617 4.86 27.3 0.74 
Conversion complete 125 39.8 0.0 —0.3014 0.0 20.1 0.66 


—90.0 








Interpretation of Results for Example Rotor 


A summary of results obtained in a complete analy- 
sis of the example rotor operated to absorb constant 
shaft power at a typical helicopter tip speed, including 
comparative data for helicopter operation, is given in 
Table 2. It may be seen that longitudinal flapping 
coefficient a; decreases steadily during conversion, while 
required 4 increases, as would be expected with de- 
creasing transverse and increasing normal flow com- 
ponents. Rotor thrust is given in per cent of thrust for 
helicopter operation. In order to indicate how much 
of the power input to the rotor is expended usefully 
during the conversion process, a profile efficiency has 
been defined as 


n = (T/P) lv — Vsin (a + a)] (6) 
or approximately, since a; is a small angle, 
n = —(T/P) (’\ + uw a)wR (6a) 


This expression represents, to a close approximation, 
the ratio of the rate of useful work done on the air 
flowing through the rotor tip path plane to the shaft 
power input. It reduces to the familiar rotor figure of 
merit in the case of hovering flight. It will be noted 
that, throughout nearly the entire conversion range of 
a, profile efficiency exceeds that for hovering. This 
indicates that there should be no detrimental effect on 
performance during the conversion process insofar as 
rotor characteristics are concerned. In general, the 
analysis indicates that a conversion process of the type 
considered will introduce no undesirable rotor behavior 
or piloting requirements and should therefore be en- 
tirely practical. 


Rotor AERODYNAMIC CHARACTERISTICS IN PROPELLER 
OPERATION 


General Considerations 


In order to obtain high-forward speeds as an airplane 
after conversion from lifting rotor to propeller opera- 
tion, reasonably high-rotor propulsive efficiency is re- 
quired. It can be shown that operating a typical 
helicopter rotor as a propeller, with no change in rota- 
tional speed or rotor geometry (except for increasing 
blade pitch to maintain constant power absorption), 
results in poor propulsive efficiency at air speeds above 


normal helicopter maximum level flight speeds. This 
poor performance is a consequence of the low disc load- 
ings that are required in helicopter design to obtain 
satisfactory low-speed and power-off flight performance. 
The condition may be explained qualitatively as fol- 
lows: 

The thrust produced by a rotor operating with con- 
stant power input is, roughly speaking, inversely pro- 
portional to the airflow velocity component normal to 
the rotor tip-path plane. For a typical rotor operat- 
ing as a propeller at an air speed of, say, 170 m.p.h., 
this normal velocity component is about ten times 
greater than in hovering as a helicopter. Rotor thrust 
is thus of the order of one-tenth of hovering thrust. 
In hovering, a representative blade section might 
operate at c, = 0.5 with a section lift/drag value of 40; 
accordingly, in high-speed flight the same section would 
operate at about c, = 0.05 with a lift/drag value of 
about 5. Had such a rotor been required to operate 
only in the high-forward speed propeller condition, it 
could have been so designed that a mean blade section 
would operate near best lift/drag with a consequent 
large improvement in propulsive efficiency at high 
speed. 

From this explanation it is clear that to improve pro- 
pulsive efficiency at high speed it is at least necessary to 
increase the mean blade-section lift coefficient in this 
operating condition. This is equivalent to increasing 
rotor power coefficient, which might be done in the 
following ways: 

(a) Increase the power input to the rotor. This is 
impractical because of the already low power loadings 
required for adequate performance in helicopter opera- 
tion. 

(b) Use contractable blades to decrease rotor diame- 
ter and tip speed during propeller operation. This 
would be an excellent solution if structural difficulties 
could be overcome. It has the disadvantage that, 
in the case of twisted blades, telescoping the blades 
would decrease overall blade twist in an operating 
régime where increased twist is desirable. 

(c) Decrease rotational speed for propeller operation. 
This appears to be the only practical solution at present. 
It requires the use of a variable-speed transmission 
and requires provision for the increased rotor torque 
realized when rotor speed is decreased. 
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Fic. 5. Envelope propulsive efficiency curves from — 


Electric Company data. a = —90°; D = 2ft.; b= 


An associated problem is provision of proper blade 
twist. In helicopter operation, moderate amounts of 
blade twist have the advantageous effects of reducing 
hovering induced power losses, reducing peak blade 
bending moments, and increasing the forward speed at 
which stall of the retreating blade tip begins to occur. 
However, it is generally agreed that the region of dimin- 
ishing returns is reached when more than about 10° 
total twist from blade root to tip is used. Although 
never demonstrated experimentally, it has been pre- 
sumed that excessive blade twist would lead to reduced 
performance in autorotation due to stalling of the 
blade root sections and, possibly, to an undesirable blade 
bending moment distribution. 

Optimum total twist for conventional propeller 
blades is of the order of 45°. The amount of twist used 
in a particular case depends on the design advance 
ratio for the propeller: the higher the advance ratio, 
the less twist required. When a propeller operates 
with nonoptimum twist distribution, propulsive ef- 
ficiency suffers because of increased profile and induced 
losses. 

The conflicting blade twist requirements for lifting 
rotor and propeller operation may be met in two ways: 
(a) use blades whose twist may be varied, or (b) 
use a compromise permanent twist distribution. The 
first alternative would be preferable if a reliable method 
for varying blade twist were available. At present, 
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Fic. 6. Comparison of calculated and experimental untwisted- 
blade propeller characteristics. 


it is thought necessary to adopt the second alternative. 
Accordingly, for the calculations of this paper unta- 
pered blades with 20° linear twist from root to tip have 
been assumed. It may be mentioned in passing that 
preliminary analyses for such blades in lifting rotor 
operation have indicated that, compared to equivalent 
untwisted blades, blade bending moments at high- 
forward speed are somewhat less severe, and that the 
initial point of stall on the retreating blade is moved in- 
board from the tip to about x = 0.65 and occurs at a 
higher forward speed. There is no evidence of impair- 
ment of autorotational performance, but the calcula- 
tions do not account for losses due to blade stall. 


Experimental Data 


Available test data on rotors operating as propellers 
are extremely limited, and, except for a few cases, pro- 
peller test data cannot be used because of wide differ- 
ences between conventional propellers and the proposed 
type of rotor in blade twist, plan form, and solidity. 


The General Electric Company has kindly made 
available unpublished results of extensive model tests 
made under its auspices. A summary of some of these 
data is presented in Fig. 5, where an envelope propul- 
sive efficiency curve obtained with a 2-ft. diameter, 
two-bladed model rotor with untwisted blades operated 
through a large range of pitch angles is shown. Pitch 
angles and power coefficients are spotted along this 
curve to indicate the operating condition. A large 
spinner with diameter approximately 25 per cent of 
rotor diameter was used on this rotor, and spinner tare 
drag was added to observed thrust in calculating propul- 
sive efficiency; therefore, resulting efficiency values are 
probably somewhat too high, since the inner portions 
of such untwisted blades normally would be stalled 
negatively in propeller operation, except at the highest 
pitch angles. For comparison, a 2-ft. diameter two- 
bladed model of a conventional propeller of uniform 
design pitch was tested concurrently, and the envelope 
propulsive efficiency curve for this model is also shown 
in Fig. 5. It should be noted that these tests were 
carried through much higher values of the advance 
ratio, V/nD, than is normal for propeller operation. 
For example, in cruising flight the DC-3 airplane oper- 
ates at about V/nD = 1.2. The conventional propeller 
model shows decided superiority in the low-advance 
ratio range, while the untwisted-blade rotor model is 
superior above V/nD = 3.0. It is notable that a pro- 
pulsive efficiency of 0.70 was obtained with the un- 
twisted-blade rotor at 80° pitch, and, for a considerable 
pitch range below this, somewhat better efficiency is 
realized. (The low Reynolds Number of the model 
tests compared to full scale probably resulted in pro- 
pulsive efficiency values lower than would be obtained 
in full-scale tests.) 


A rapid increase in power coefficient required for 
operation at envelope efficiency as pitch is increased is 
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Fic. 7. Comparison of actual and assumed equivalent blade 
pitch distributions. 


apparent. As a consequence, operation with constant 
power input would require that tip speed for % = 80° 
be only 40 per cent of the tip speed required at % = 
35°. These results indicate that even an untwisted- 
blade rotor might perform satisfactorily as a propeller 
under proper operating conditions. 

In Fig. 6 test results for the untwisted blade rotor at a 
particular pitch setting are compared with values cal- 
culated by the method described below. The fair agree- 
ment indicates that the calculation method provides a 
reasonably close prediction of magnitude and trend. 
Part of the discrepancy may be explained by the fact 
that the rotor assumed in the calculations had lower 
solidity and blade sections operating at a much higher 
mean Reynolds Number than the model rotor, which 
would tend to reduce induced and profile losses for the 
calculations in comparison to the test data. 

A second source of applicable data is found in refer- 
ence 1, where excellent data are presented from tests of 
34-in. diameter model propellers, including several 
having unusually low amounts of blade twist and 
roughly constant chord from shank to tip. The twist 
distributions for two of the model blades are shown in 
Fig. 7; for purposes of this analysis it has been assumed 
that test results for propellers using these blades would 
be closely approximated by untapered blades having an 
equivalent linear twist distribution given by the dashed 
lines of the figure. From Fig. 1 of reference 1, the 
estimated effective value for the solidity of these propel- 
lers is ¢ = 0.120. This is about triple the solidity 
desired for the type of rotor under discussion. How- 
ever, it has been assumed that, for a given operating 


condition, thrust and torque coefficients vary ditectly 
with solidity, and propulsive efficiency is independent 
of solidity. These assumptions are valid only for 
conditions where axial and rotational induced velocities 
are negligibly low, which conditions are approximately 
fulfilled in the present cases. 

In Fig. 8 propulsive efficiency curves are shown for a 
typical low-solidity rotor operating at several different 


‘power coefficients. These curves were interpolated 


from the test data of reference 1 in accordance with the 
assumptions stated above. In order to show their 
physical significance, the Fig. 8 curves have been re- 
plotted in Fig. 9 in terms of the operation at various 
tip speeds of a 19-ft. diameter rotor absorbing 100 
shaft hp. The curve for 671 ft. per sec. tip speed 
shows propeller operation at a tip speed that might be 
used for a typical modern rotor in helicopter operation, 
Propeller performance is poor in this condition. When 
tip speed is reduced, a considerable increase occurs both 
in peak efficiency and in air speed for peak efficiency. 
The curve for 230 ft. per sec. tip speed represents accept- 
able propeller performance for a moderately high-speed 
airplane. 

In Figs. 10, 11, and 12, some of the test curves of 
reference 1 for constant pitch settings are shown as 
interpolated for 20° linear blade twist. (The nominal 
pitch settings of reference 1 have been increased by 3° 
to account for the mean difference between nominal 
and zero-lift blade chord lines as indicated by Fig. 3 of 
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Fic. 8. Propulsive efficiency curves for operation at constant 
power coefficient, based on model test data. o = 0.04; 4 = 
—20°; a = —90°, 
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Fic. 9. Effect of tip speed on propulsive efficiency for con- 

stant power operation, based on model test data. a = —90°; 

D = 19 ft.; « = 0.04; 0 = —20°; b.-hp. = 100; altitude = 
sea level. 
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Fic. 10. Comparison of calculated and model test thrust coef- 
ficient variation with advance ratio. a = —90°. Test 
data, 34 in. diameter, o ,0.12, interpolated for 6, = —20°. 
—o— Calculated, o = 0.04, 6, = —20°. 
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Fic. 11. Comparison of calculated and model test power coef- 
ficient variation with advance ratio. a = —90°. -~——-Test 
data, 34 in. diameter, o = 0.12, interpolated for 0, = —20°. 
—o— Calculated, « = 0.04, 0, = —20° 
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Fic. 12. Comparison of calculated and model test propulsive 
efficiency variation with advance ratio. a = —90° Test 
data, 34 in. diameter, o = 0.12, interpolated for 6; = —20°. 
—o— Calculated, « = 0.04, 0; = —20°. 


reference 1.) These are compared with values calcu- 
lated by the method described in the next section. 
Agreement is fair, although differences appear to in- 
crease with increasing pitch. The calculations yield 
peak propulsive efficiencies about 10 per cent too high. 
As in the case of the General Electric Company data, 
part of the discrepancies may be explained by the fact 
that the rotor assumed in the calculations had a much 
lower solidity and had blade sections operating at a 
higher mean Reynolds Number than the model propel- 
ler. Also, hub drag was not accounted for in the cal- 
culations. On the basis of these comparisons, it is 
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probably correct to conclude that the calculation 
method will give reasonably accurate results for the 
somewhat unusual operating conditions required for 
the rotor-propeller. 


Calculated Aerodynamic Characteristics 


The method of propeller analysis used for the calcula- 
tions of this paper is an adaptation of the methods of 
references 2 and 3, which are based on standard vortex 
theory with Goldstein corrections. In reference 3 a 
comparison between experiment and calculations by 
this method is shown, in which excellent agreement was 
obtained through V/nD = 6.0 in spite of blade load- 
ings that varied widely from that on which the Gold- 
stein theory is based. In the present case, blade load- 
ings vary more widely from that of the Goldstein theory 
than did the loadings of reference ‘3; however, it is 
thought that reasonable agreement should still be ob- 
tained. 

A diagram of the flow condition at a typical blade 
element is shown in Fig. 13. That part of the induced 
velocity which is due to profile drag is neglected. €isa 
small angle. 

The fundamental relation of the vortex propeller 
theory as modified by Goldstein may be written as 


= BcWo,/S8nrrF sin o (7) 
or, 
tane = w/W = oc¢,/8xF sin > (7a) 


From the diagram, Fig. 13, 





W = (wr/cos ¢) — w tan ¢ 
= (wr/cos ¢) W(tan e) (tan ¢) (8) 
W sin @ = wr/(cot @ + tan e) (9) 
- 
Wrur 
b 

















Fic. 13. Schematic diagram of a ——. blade section in 
propeller operation. a = 
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Then, 


dCr - ] qr — . tae x 
dx (p/2) (wR)*eR? dr ~~ (p/2) (wR)? rR? 
(p/2)Wbc(c, cos 6 — Ca sin >) 


cot @ — tan k 


= S§Fx? (tan e) = (10) 
(cot @ + tan e)? 
dCp _ ] dP - a 
dx (p/2) (wR)* eR? dr (p/2) (wR)? rR? 
(p/2) Wbcr (c¢, sin @ + Ca cos >) 
= 8Fx‘ (tan e) Lt 28 te * (11) 





(cot @ + tan e)? 


Using the above expressions, calculations were made 
for a range of blade radial stations and pitch angles, 
and a range of section angles of attack at each combi- 
nation of radial station and pitch angle. Aerodynamic 
characteristics for the NACA 23018 blade section esti- 
mated from data of references 4, 5, and 6 are shown in 
Fig. 14. The accuracy of the method is not sufficient 
to justify consideration of Reynolds Number variation 
with section radius. A sample calculation is shown in 
Table 3. Results for each radial station were plotted 
as curves of C,/¢ and C,/o versus V/nD for constant 0. 


Having a complete set of blade-element characteris-’ 


tics curves for blade stations x = 0.15, 0.30, 0.45, 
0.60, 0.75, 0.85, and 0.95 and for the operating range of 
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Fic. 14. Blade section characteristics used for propeller per- 
formance analysis. 
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Typical calculated thrust coefficient grading curves. 
Oz = 0.7%) = 39°; 8 = —20°; o = 0.04. 
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TABLE 3 

Sample Calculation of Blade Element Characteristics 
(i) = 0.60 
(2) 6 (deg.) 45.0 

(3) a» (deg.) (assumed) 8.0 

(4) c (Fig. 14) 0.77 
(5) tan Rk (Fig. 14) 0.0230 
(6) @¢ = 6 — a (deg.) = (2) — (3) 37.0 

(7) sin @ = sin (6) 0.602 
(8) F(6 = 3) (from reference 2) 0.770 
(9) o 0.040 
(10) tan e = o¢:/8xF sin @ = (9) (4)/8(1) (7) 0.01384 
(11) ¢€ (deg.) = 57.3 (10) 0.79 
(12) oo = @ — e (deg.) = (6) — (11) 36.21 
(13) V/nD = xx(tan ¢@) = m (1) tan(12) 1.381 
(14) cot @ = cot (6) 1.327 
(15) cot @ — tank = (14) — (5) 1.304 
(16) (cot @ + tan e)? = [(14) + (10)]? 1.800 
(17) 8Fx* tan e = 8(8) (1)% (10) 0.01840 
(18) dCr/dx = (17) (15)/(16) 0.01333 
(19) 1 + (cot ¢) (tan k) = 1 + (14) (5) 1.0305 
(20) 8Fx‘ tan e = 8(8) (1)4 (10) 0.01105 
21) dCp/dx = (20) (19)/(16) 0.00642 


section pitch angles, thrust and torque coefficient grad- 
ing curves were constructed for the desired pitch set- 
tings and blade twist distributions at intervals of ad- 
vance ratio to provide comparisons with the model test 
data. A typical set of these grading curves is shown in 
Fig. 15. The areas under the grading curves corre- 
spond to the thrust or torque coefficients for their re- 
spective operating conditions. 

Comparisons of calculated and test results are shown 
in Figs. 6, 10, 11, and 12. These have been discussed 
previously. 


EFFECT OF GRAVITY ON A LAG-HINGED ROTOR IN 
PROPELLER OPERATION 


Introduction 

In the case of a rotor with lag hinges operating as a 
propeller with rotor axis horizontal, the alternating 
gravity force on each blade may be expected to produce 
a blade motion response. The following analysis in- 
vestigates this blade motion and the resulting forces 
on the hub. 


Assumptions 

The following simplifying assumptions are made: 

(a) y is a small angle, and y? terms are negligible 
compared to y terms. 

(b) For purposes of this analysis, blade aerodynamic 
forces may be neglected, since, because of symmetry of 
the operating condition, such forces are essentially con- 
stant. 

(c) The rotor hub rotates at constant angular veloc- 
ity. 

(d) The effect of lag motion damping may be neg- 
lected in considering the limiting case. 


Blade Lag Motion 


Fig. 16 shows a front view of one blade and rotor 
hub in propeller operation. The blade lag motion is 
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analyzed with respect to coordinates rotating with the 
hub. Centripetal acceleration is then replaced by a 
centrifugal force given by 


F, = (W,/g) w*R(e2 + Li) (12) 
plus Coriolis terms containing y which become negli- 
gible in Eq. (13) since they are of y* order. The equa- 
tion for angular motion about the lag hinge is 


2y yeoR 
d Y aes F.LoR ye2 


—_ = — WLR si ~ = 
dt? eR + LR LR sin (y 7) 


I: 


wR? 2 L 
—-W, | Pe Pin Fe th. y¥ +L, R(sin Y — y cos »| 
g eo + Le 


(13) 
Assume that the blade lag motion is given by 
y¥=A,cosy + B,siny + Az cos 2y + Bzsin2y (14) 
Substituting Eq. (14) in Eq. (13) yields 
Qj = KA, - COS y (KiA, ad KA») + sin y [KiB “+ 
K2(Be _ 2)] oe cos 2y (KiA2 + 312A» + K2A}) a 
sin 2¥(K,B. + 3/2B. + K2B,) + cos 3y (K2A2 +...) + 
sin 3¥(A.B. +...) (15) 
Equating coefficients of the constant, sin y, cos y, and 
sin 2) terms to zero (since there are only four unknowns, 
only four of the seven possible simultaneous equations 
are solved, and these are selected under the assump- 
tion that the higher harmonics will prove to be rela- 
tively unimportant), 
A, = 0 (since Kz is not, in general, zero) (16) 


Az = 0 (since K is not, in general, zero) (17) 


B, = (K2/K;,) (2 — Bs) (18) 

B, = —B.(K, + 3l2) Ke (19) 
Solving for Bo, 

B, = 2/[1 — K,(Ki + 3k2)/K:?] (20) 


Motion of Blade Center of Gravity 

Referring to Fig. 13, the motion of the blade cen- 
ter of gravity in fixed coordinates is derived as fol- 
lows: 

y = eR sin y + LR sin (yy — vy) 

~ R[e. sin y + L, (sin yy — ycosy)] (21) 
z= eRcosy + LR cos (vy — y) 

~Rle cosy + Li (cosy + ysiny) (22) 


Substituting y = B, sin y + B, sin 2y and differentiat- 
ing twice, 


d*y 
] 2 


adi- 


= w*R}{sin y[ — e — Li + ('/2) LiBe] + 


sin 2¥(211B;) + sin 3y [(9/2) L:Bs]} (23) 
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Fic. 16. Schematic diagram of lag-hinged rotor blade in pro- 
peller operation. 
d*z 


it? - w*R } cos v[— e — Ly — ('/2)L, Be] + 
at 


cos 2p [2L,B,] + cos 3p[(9/2)L:Bs]} (24 


Hub Loads Due to Gravity 

The equations of linear motion of the blade center of 
gravity are 

- W, = (W,/g)(d’2/dt?) (25) 
F,, = (W,/g)(d*y/dt?) (26) 
where the vertical force exerted by one blade on the hub 


at the lag hinge is F,, and the horizontal force is F,,. 
Substituting Eqs. (23) and (24). : 


F,, = W, (1 + (w*R/g) {cos y[— ae — Li — 

(1/2)L,Be] + cos 2¥(2L,B;) + cos 3p [(9/2)L1Bs] } ) (25a) 

F,, = (Wyw?R/g) {sin y[—e2 — Li + (1/2)L:Be] + 
sin 2¥(2Z,B,) + sin 3p[(9/2)L,B.]} (26a) 


For a three-bladed rotor, the total instantaneous blade 
reactions on the hub are 


FP, = F,, + F,, + F,, 27) 
Fi = F, + F,, + F,, (28) 


where the value of ¥ for each of the three blades may 
be expressed in terms of y for the first blade as 


W=v; ve=Wvt+ (27/3); vs = w+ (42/3) 
Making these substitutions, 
27 Www? R en 
a ————— Best +3W, (0) 
g 


Me 5 
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; 27 Wywe?L.R 
F, = —— 
2 g 





B, sin 3y (30) 


By an analogous process for the case of a two-bladed 
rotor, 


4Wyw? LR 
ji, = ————— Bre 9 + 2, (31) 
4Wyw*LR 
Fy = — — B, sin 2p (32) 
g 


In the case of a three-bladed rotor, the force transmit- 
ted to the hub by the blades in propeller operation 
(aside from a steady vertical force equal to total blade 
weight) is a third harmonic vibratory force arising 
In the case of a 
two-bladed rotor, the transmitted force is a second 


from the second harmonic lag motion. 


harmonic vibratory force arising from the first harmonic 
lag motion. 


Numerical Examples for Typical Rotors 


The foregoing analysis has been applied to four 
typical rotors, selected to represent extremes of size 
and operating at low rotation speed (found to be more 
critical in the case of three-bladed rotors) to represent 
propeller operation. Results are summarized in Table 
4, In the case of the three-bladed rotors, the vibra- 
tory hub loads are negligibly low compared to total 
blade weight, while in the case of the two-bladed rotors 
the vibratory hub loads are dangerously high. These 
results indicate that a three-bladed rotor with lag 
hinges should be satisfactory from the standpoint of 
gravity-induced vibratory hub loads; however, for a 
two-bladed rotor, a semirigid or rigid-type rotor would 
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SCIENCES—OCTOBER, 1949 
TABLE 4 = 
Rotor diameter (ft.) 16.5 16.5 72.0 72 | 
Number of blades 3 2 3 2 
Weight of one blade (lbs.) 10 15 200 300 
Rotor speed (rad. per sec.) 32 32 8.2 82 
B, (rad. times 107%) 7.75 7.78 ST .7 57.7 
By (rad. times 10~) —6.50 -—6.50 —223.0 —293 
Vibratory hub load result- ; 
ant amplitude (Ibs.) 0.15 44.2 23.9 2745 





be required to avoid objectionable vibratory loads in. 
duced by gravity. 
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Matrix Solution of the n-Section Column 






WILLIAM T. THOMSON* 


University of Wisconsin 


SUMMARY 


A procedure is developed for the determination of the critical 
load of a column of » different sections under any end conditions. 
When 7 is a small number, the method leads to simple analytic 
expressions for the critical load 


NOMENCLATURE 


l = length of ith section 

E;l = stiffness of 7th section 

Xx = x coordinate at the end of the 7th section 
yi = deflection at the end of the 7th section 
yi’ = slope at the end of the 7th section 

3 = V P/E: 

E = Young’s modulus of elasticity 

P = column load 

l = total length of column 

Mo, M, moment at ends of the column 


DEVELOPMENT OF EQUATIONS 


Fr: | SHOWS A COLUMN of » sections with arbitrary 
end conditions. 


so that failure tends to take place by lateral instability 


The column is assumed to be long 


rather than by direct compression. 


From the free-body diagram of the ith section shown 
in Fig. 2, the differential equation for this section be- 


comes: 


El ,(d*y/dx*) = —P(y—y-1) + Mi-i (D 


Employing the method of Laplace transformation, the 
solution for the deflection at 7 can be written in terms 
of quantities at 7 — 1 as 

sin Bl; M,;-1 


Vi =M-1tM-1 (1 


: P cos Bl; ) 
B; Ed iB; 


(2) 


Differentiating, the slope and moment equations are ob- 


tained: 
, , M; =e ‘ 
yi = ¥i:-1' cos Bd; + = sin Bl; (3) 
? “it iPt 
M, = —y,-1'Ed iB; sin Bd; + My—10cos Bd; (4) 


Expressed in matrix form, Eqs. .(2), (3), ‘and (4) are 
represented as follows: 





y 1 (1/8,) sin Bd; (1/EWd,8)(1 — cos B,l,) Sims 
y;,/ |} =10 cos Bd; (1/E,J;8;) sin Bd, Vi'-1 (5) 
M, 0 EJB, sin Bil; cos Bl; M; -1 
Itis evident then that by repeated multiplication of the (2) Fixed-end columns: 
square matrix the quantities at m can be expressed in 7 See soe 
terms of corresponding quantities at 0 by an equation of eet te SF 
the form: Therefore, 
Aj; = Ax = 0 (8) 
Yn Ai Ay Aj3 Yo ; : ; , 
vy,’ | = 0 as Me Yo! (6) (3) End 0 fixed and end n pinned and free: 
M,, 0 Ay Ass Mo yo= y = M,=0 
Knowing the end conditions, the equation for the criti- Therefore, A 
3s = 0 (9) 


cal load is readily determined from Eq. (6) as fol- 
lows: 


(1) Pin-ended columns: 


Yn => Vo = M, = M, = 0 


Therefc re, 


Ax = Azo = 0 ( 


~I 
~— 


Received April 4, 1949. 
* Associate Professor of Mechanics. 


In interpreting Eqs. (7) to (9), when two of the ma- 
trix elements are equated to zero, only roots satisfying 
both equations represent the true solution. 


PROCEDURE FOR COMPUTATION 


The matrix elements of Eq. (6) are obtained by mul- 
tiplying the m square matrices of Eq. (5). Whenmnisa 


large number, it is simpler to carry out this computation 
numerically as follows: Choose a value for 8,, and de- 
termine 
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Mo 


Fic. 1 (left) and Fic. 2 (right). 


, B,-2= Bn; a 


etc. Compute the elements of the square matrix of 
Eq. (5) for each section, and multiply out to obtain 
the elements of Eq. (6). Plotting Eqs. (7), (8), or (9) 
as function of 8,, the critical load is determined from 
the particular values of 8, where the curve goes to 
zero. 

When 1 is a small number, simple analytical expres- 
sions can be obtained for the critical load. Some of 
these are given in the following section. The case for 
n = 1 is included since these equations are well known. 


ANALYTICAL EXPRESSIONS FOR SMALL 1 


The end conditions are designated as (1), (2), and (3) 
to conform to Eqs. (7), (8), and (9). 
Forn = 1.0: Case (1), 


sin fi = 0, Bl = VP/Eli=-¢ 


SCIENCES—OCTOBER, 1949 


Therefore, 
P.. = wr°EI /|? 10 


Case (2), 


(1 — cos @l/) = sin pl = 0, Bl = 2x 
Therefore, 
Pi, = 4n°EI /I? 1] 
Case (3), 
cos pl = 0, pl = r/2 
Therefore, 
P., = wEI/Al? 12 
Forn = 2: Case (1), 
_ ee = ro a aah) ai 
sin Bols sin §;/, 


- ; cos B,/; — COS Bol. = 0 


Bi 
Therefore, 
tan Bolo = —(82/B1) tan By), (13 

Case (2), 
1 — cos Bl; cos Belo + (81/ Be) sin Byl; sin Bolo = Ol 14) 
8, COS Bole sin 8y/; + Be sin Bole cos Bil, = O ho 
Case (3), 

— (B,/B2) sin B,J; sin Bolo + cos Bil; cos Bole = 0(15 

Forn = 3: (Case (1), 

tan 8); 

Wea 


—_ (—— —- 16 


Bo 
; — — tan §,/; tan aul | 
8; 


tan Bol | 
tan B31 
a 





Case (3), 
— (8,/82) tan Bol. tan 6]; + 1 = 0 (17 


The correctness of these equations can be easily 
verified by letting 8, = 62 = 6; = ete. =) =); 
= ]/n in which case the results reduce to that of n = 1. 


CONCLUSION 


A matrix method is developed by which the critical 
load of an m-section column can be determined for any 
end conditions. For m a small number the results 
reduce to simple equations for the critical load. 
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A Modification of the Holzer Method for 
Computing Uncoupled Torsion and 





Bending Modes 


HENRY E. FETTIS* 
dir Materiel Command 


SUMMARY 


The present method is an adaptation of the well-known process 
devised by Holzer for solving the differential equations of beam 
vibrations. Several versions of this method are now in common 
use; the original form-is described in references 1 and 2 and is 
usually employed in computing the torsional modes of vibration 
of a bar carrying a series of concentrated discs. Similar processes 
were devised by Myklestad® and Bellin® for bending modes, and 
recently an extension of the above principle to the computation 
of coupled modes has been published.4 All of the above methods 
have in common the feature that a value of the frequency is as- 
sumed and a deflection computed which satisfies all but one of the 
Although the same principle forms the 


basis for the following method, the system of numerical integra- 


required end conditions. 


tion employed makes the method better suited to configurations 
with continuous mass and stiffness distributions. The process 
is set up in such a manner as to be easily adapted to routine 


calculation 


by a semiskilled computer and the correctness of 
the results may be verified at each station. All of the features 
that make the Holzer method practical are retained, and some 


new techniques not in common practice are introduced. 


1) FORMULATION OF THE METHOD 
(1) Torsional Modes of Vibration 


p | {iE DIFFERENTIAL EQUATION to be solved is: 


d/dx)GJ(d0/dx) = —Iw*6 (1) 


where @ is the angular deflection in radians in a normal 
mode of vibration, w is the frequency in radians per 
second, GJ is the torsional stiffness in pound feet’, J is 
the mass moment of inertia per unit span about the 
torsion axis in pound seconds’, and x is the spanwise 
coordinate in feet measured positive to the right along 
the torsion axis. 

Eq. (1 
simultaneous equations: 


*Xn.1 
GJ, +. 100 +. 1 a GIB,’ _ —a f 6 dx 
Hn (2) 


Xn +1 
8 Ao 6, -f 6'dx 
Xn 


in which subscripts indicate the value of the variables 
at any two points x, and x, 4 ,0n the span (x, < X, +1) 
and primes denote differentiation with respect to x. If 
the interval 


may be written in integrated form as two 


AX, = Xu 4.1 — Xe 


Received March 11, 1949. 
* Mathematician, Dynamics Branch, Aircraft Laboratory. 


is small, the integrals may be approximated by the 
trapezoidal rule, and Eqs. (2) become 


: Ax, 
1 — GJ,’ = —o( : x 


(i, + 10, + 1 + T An) (3a 


; AX», , : ; 
6.4+,—- 0, = : (0, 41° + 0,’) (3b) 


GJ, + 19, 


6, + 1 may now be eliminated between Eqs. (3a) and 
(3b), giving a simple linear relation between the de- 
flection at any point (or “‘station’’) and the deflection 


and slope at a preceding point: 


: ,f Axa\? 
On +1) GJ, 41 + -— In 41 - 
Ax, 2 
0] Cl 2. = u*( Pa ) | + (GJ, as a GJ) x 
»[{ AXn 
6, - (4a) 


The slope at the advance point may then be found by 
rewriting Eq. (3b) in the form: 


On 41° = (2/ Axn) (On +1 — On) — 8,’ (4b) 


Eqs. (4a) and (4b) may be employed to compute a 
normal mode of vibration according to the usual Holzer 
scheme by assuming a value for the frequency w and 
starting the solution at one end of the beam where 
either the deflection or slope or a linear relation between 
them is known. These initial values are determined 
by the end conditions; from physical considerations, 
the following end conditions will apply: 
At a free end carrying no concentrated mass, 


GJ(6') = 0 (.) 
At a fixed end, 
¢é=0 (6 
At a free end carrying a concentrated item of mass 
moment of inertia J, 


GJ(0’) = =wl@ (7 

+ The bar over the letter J is here employed to distinguish be- 
tween a concentrated moment of inertia and a distributed inertia 
as in Eq. (1 
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and sign being plus or minus according to whether the 
end is the right or left one.* 


(2) Bending Modes 
In this case the differential equation is: 

(d?/dx?) EI(d*h/dx*?) = pw*h (8) 
in which A is the deflection in feet, 4 is the mass per 
unit length in pound feet~? seconds’, EJ is the bending 
stiffness in pound feet’, and w is the frequency in sec- 


onds~!. Eq. (8) may be written in integrated form as 
four simultaneous equations: 


Xn +1 
| Fe4i1-F, = wf ph dx (9a) 
Xn 

| i 

ae me I F dx (9b) 

4 Xn 

| *Xy +1 

| Pn +1 ~~ Pn _ J (M/EI1 \dx (9c) 
Xn 

| X,, +1 

Reaai-k = f p dx (9d) 
Xi 


which, when the integrals are approximated by the 
trapezoidal rule and the quantities F, 41, \/, 41, and 
bn +1 eliminated, give the following relation between 
the deflection, h, 4 1, at any station and the deflection 
h,, slope p,, shear F,, and bending moment .\/, at the 
preceding station: 


: »{ d*n\' 
hy, + {1, oe eae (5 ) Kn + 1 = 
‘ of Axn\' Az,\.. 
hy| El +1 + o*( ‘a ) a + 2a > et, ‘2 + 


M,f Ax,\?._. ; — ( Axn\3 
EI ( ~ ) [EI, +1 + El,] + 2F,( = ) (10a) 


Advance values of p, M, and F are then found by re- 
arranging Eqs. (9b), (9c), (9d): 


ax, ax. " 
| Paty ) = [hn +1 — ha] — Pal = (10b) 


| M,, +1 Ax, “ 1 Att, MM, 
| . ¢ [Pr + 1 r Pnj 7 = : 
EI, 41\ 2 2 fils 


| + 


. AX» ° AX , 

| a , ( ) = [M,, 5 Mall > ) 

2 , 2 

| (Ax, \3 

Fa = ) (10d) 


From these equations, the four quantities /, p, F, and 
M may be computed when their values at a preceding 
station are known. To start the solution, use is made 
of the boundary conditions at one end appropriate to 
the physical conditions of the problem. End conditions 
applying in conventional cases are as follows: 


(10c) 


ll 


* This sign convention holds as long as x is measured positive 
to the right. 
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At a free end carrying no concentrated mass, 


Elk” = 0 (11a) 
(EIR")’ = 0 (11b) 
At a clamped (built-in) end, 
fh=0 (12a 
lh’ =0 (12b 


At a simply supported (hinged) end, with no concen- 
trated moment, 
(13a 
(13b 


jh =0 
\(EIh") = 0 
At a free end carrying a concentrated item of mass m 
and mass moment of inertia about a chordwise axis 
(rotary or bending inertia) equal to /, 
f +(EIh") = —w'Ih’ 
( +(EIh")' = w*mh 


(14a) 
(14b 


the sign being plus or minus according to whether the 
end is the left or right one.t 
At a hinged end with a concentrated item of rotary 
inertia J, 
frh=0 (15a 


)=(EIh") = —wtTh' 15b 


The problem in the bending case is complicated by 
the fact that two of the four variables must be made 
arbitrary at the start of the process. Thus, each suc- 
cessive value of the deflection slope, moment, and shear 
must be written as a linear combination of the two 
arbitrary initial values. 

Details of the two above processes are worked out 
and exemplified in Section (III). 

(II) APPLICATION TO AIRCRAFT WINGS: ADVANTAGE 
OVER OTHER COMMONLY USED METHODS 


The modes of vibration of aircraft wings may evi- 
dently be approximated by the modes of a free-free 
beam carrying a concentrated load at the center. It 1s 
simpler, however, to consider only the semiwing from 
the centerline outward and to divide the problem into 
the two cases of vibrations that are symmetric with 
respect to the aircraft plane of symmetry and _ vibra- 
tions that are symmetric with respect to the centerline. 
The first type are commonly referred to as ‘‘symmetric 
and the latter as ‘‘antisymmetric’’ modes of vibration. 

The end conditions summarized in Section (1) are 
applicable to the present problem with only slight modi- 
fications. Thus, at the wing tip, conditions (5) or (7) 
apply for torsional modes, while Eq. (11) or Eq. (14) is 
used in the bending case. At the airplane centerline, 
condition (6) holds if the torsion is antisymmetric, 
and condition (7) holds if the mode is symmetric tor- 
sion. In the latter case, / is to be interpreted as half 
the moment of inertia of the fuselage in pitch. For 
antisymmetric bending, the appropriate conditions are 


t See footnote * on this page. 
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COMPUTATION OF TORSION 


(15) with] replaced by half the moment inertia of the 
fuselage in roll, while for symmetric modes the condi- 


tions become 
k’ = 0 (16a) 
aeury = —mwh (16b) 
where m is half the mass of the fuselage. 

One of the principal advantages of the Holzer method 
lies in the fact that higher modes may be computed 
without reference to the lower ones and that, by choos- 
ing the point x = 0 at the wing tip, both the symmetric 
and antisymmetric cases may be handled simultane- 
ously, since the two cases are distinguished only by the 
boundary conditions at the center. Other methods 
require separate treatment in the two instances, and 
little or none of the information obtained from one 
can be used in the other. 

An additional advantage gained by the present 
treatment of the problem lies in the fact that discon- 
tinuities in the basic parameters which frequently oc- 
cur, due to structural changes, cutouts, or large concen- 
trated masses at inboard points, are more easily ac- 
counted for provided the stations are so selected that 
they coincide with any discontinuities that may exist. 
Consider, for example, the case where the torsional 
rigidity changes suddenly at some point, x,. -At the 
point, GJ is two-valued, having the value (G/J;,), when 
approached from the left or outboard side and the value 
(GJ,)z when approached from the right or inboard side. 
Since the torque across this section must remain con- 
tinuous, it follows that (@’) must also make a finite 
jump at the point. The slope on the left side is found 
from the usual equations, viz., 


0.) (=. ' = 1, — -:] — & Y oo ‘) 
(17) 
whereas the slope on the right will be determined from 
the condition of continuity of torque—namely, 
(GIy) (Oe )z = (GIe) R(x’) p (18) 


The angular displacement remains continuous at the 
point, and so the advance value of 6 is given by 


: Jd Bx,\* 
el Gh +1 + «'( ° ) I, | _ 
, , Ax; 2 
af Gu A es u*( 5 ) i | + 


Axx 
isda + Glos 100") =) (19) 


A similar situation occurs when a concentrated mass is 
located at an inboard point. In this case, the torque, 
and therefore the slope, is discontinuous and the slopes 
on either side of the point are related by the equation 


GI, [ (0.2 = (6: ) | = Tw, (20) 
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where J is the mass moment of inertia of the concen- 
trated item. Eq. (7) may be regarded as a special 
case of Eq. (20), since then either (6,’), or (®’)p is 
zero, depending on whether the end is the left or right 
one. 

In the event that the distributed inertia of the wing 
is discontinuous, a different value of J must be used 
when computing the displacement at the discontinuity 
than is used when proceeding from this point to the 
next one. However, no discontinuity in either angular 
displacement or slope will result in this case (cf. Ex- 
ample 1, Section IV). 

Analogous situations to the above in the bending 
case may be disposed of in a similar manner. Thus the 
presence of a concentrated item of mass m and rotary 
inertia J produces a discontinuity in both shear and 
bending moment, such that 


§ (Fr)z onl (Fide = —mw*ly 
war — (Me = Iw*p, 


with the slope, p, and displacement, #, remaining con- 
tinuous. Eqs. (21) include Eqs. (14) as a special case 
where either F,; and M, or Fp and Mp are zero. 

The present scheme is also useful when the effect of 
the variation of the mass or inertia of some concentrated 
item is desired. Since the mass (or inertia) of such 
items enters linearly into the equations, either by virtue 
of the end conditions or the conditions of continuity, it 
follows that the deflections, slopes, etc., at all points in- 
board of the concentrated item will be linear functions 
of the mass or inertia. Introduction of this additional 
variable into the problem permits the end conditions 
at the centerline to be satisfied for any assumed fre- 
quency. Thus, in the torsion case, the end condition 
at the center may be satisfied by proper choice of the 
inertia of some concentrated item, while in the bending 
case the two conditions at the center will result in two 
simultaneous equations involving the mass and the ar- 
bitrary tip slope. As a special case, the centerline 
condition for symmetric torsion may be satisfied by 
proper choice of the inertia at the center. Similarly, 
for bending modes, the mass or inertia at the center 
For certain 


(21) 


may be employed as an auxiliary variable. 
frequencies, negative or infinite masses may result. 
The negative values may be interpreted as effective 
masses due to flexibility [cf. Eqs. (25), (31)], while 
infinite values indicate a node at the point of attachment 
of the mass. 


(111) TABULAR FORM FOR ROUTINE COMPUTATIONS 


(1) Symmetric and Antisymmetric Torsion 


(a) Tabular Form.—See Table 1. 

(b) Outline of Procedure—Columns (7) and (10) are 
first computed from the basic data listed in columns (2) 
and (3). Columns (5) and (8) are then computed for 
an assumed value of the frequency, w. 
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TABLE 2 


Tabular Form for Bending Computations 





The initial value in column (6) is obtained from the 
boundary condition at the free end—namely, 
(GJo)(09’) = —I,w7O (22) 


where J, is the mass moment of inertia about the torsion 
axis of the concentrated item at the tip. 

The initial value in column (4) [and therefore % 
in Eq. (22)] is arbitrary and may be taken as unity. 


Succeeding values for columns (4), (9), and (6) are 
obtained from the expressions: 
(4)n41 = [An X (5)n + (On X (7), | (S)») 
(9) 3. > [(4) n “a ( +), ] _ 6)» (23) 
(6) » +1 = (9) n a x (10), + 1 \ 


continuing until the deflection and slope at the center- 
line (base) have been found. At the base, the correct- 
ness of the assumed frequency is determined from the 


condition 
(24a) 
for antisymmetric modes, and the condition 


Iyw?0, = (GJ,)(0y') (24b) 


for symmetric modes. If one or the other of these con- 
ditions is not satisfied, the process is repeated tising a 
new assumed frequency. After several trials, a plot of 
the base values, Eqs. (24a) and (24b), as a function of 
the assumed frequency indicates the location of the 
correct frequency (or frequencies). In the symmetric 
case, Eq. (24b) may be solved for J,, and the resulting 
value may be plotted as a function of the assumed fre- 
quency. The latter method enables obtain 


solutions corresponding to various loading conditions 


one to 


of the aircraft. 
The tip inertia may also be used as an auxiliary vari- 


able. In this case, the initial value of coiunin (6) is 
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written as a factor times the inertia /, of the tip item, 
Each succeeding value of @ and @’ then becomes a cop. 
stant plus a second constant times/;. At the base, the 
equation for either the symmetric or antisymmetric 
mode may be solved for I, Thus, any assumed fre. 
quency determines a unique value for the tip inertia 
in each mode, with the possible exception of those fre. 
quencies that cause the coefficient of , to vanish. Such 
a condition could be interpreted as corresponding to a 
fixed tip and would thus be of little interest in the pres. 
ent problem of aircraft wing vibration. However, in 
other cases of torsional vibration, this condition may 
also have practical significance. 

The negative inertias that will result for a certain 
range of frequencies have meaning if flexibility of the 
tip item is allowed. In this case, the effective inertia 
acting at the tip may be shown to be 


e [1 = 1 w,)"| (25 


in which w, is the natural frequency of the item on its 
restraining spring. The quantity defined by quantity 


(25) will evidently be negative if 
We < w 
that is, provided the spring is sufficiently weak. 
In a manner similar to the one described above, the 
inertia of a concentrated item located at an inboard 
station may be used as an auxiliary variable. 


Symmetric and Antisymmetric Bending 


(2) 

(a) Tabular Form.—See Table 2. 

(b) Outline of Procedure-——Columns (7), (9), (11), 
(13), and (14) are first computed from the basic data 
listed in columns (1), (2), and (3). Columns (5) and 
(12) are then computed for an assumed value of the 
frequency, w. 

The initial values in columns (8S), (10), and (15) are 
obtained by making use of the boundary conditions at 


the tip—namely, 


Mo = —I wo, Fy = mM who 


where #7, is the mass and J, is the mass moment of in- 
ertia in bending of the tip item. The initial values 
in columns (4) and (6) are arbitrary, and one of them 
Each succeed- 


(usually /y) may be set equal to unity. 
ing entry in columns (4), (6), (8), (10), and (15) thus 
consists of two parts: one depending on the initial dis- 
placement at the tip and the other on the initial tip | 
angle. These successive entries are found by making | 
use of Eqs. (10), viz., 


(An air = [Ba X (5)n + (On X Wa + (S)n X 
(9)n + (10), K (11)a)/(12) 

(6)n 41 = (14)n 41 X [An 41 — Adal — (© 

(14), ii” x [(6) io (6) n] aaa (8) n 

(15)n41= (S)n 41 X (13)n 41 

(10)n41= (14)n 41 X [(15)n 41 — (15)0] — 


(S)n 7 > 


(10), 


(26 
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COMPUTATION OF 


(n = 0, 1,2, ...). At the base, the deflection, slope, 
moment, and shear are thus obtained in terms of the 
arbitrary tip deflection and arbitrary tip slope. The 
poundary conditions at the base may then be applied 
to determine the correctness of the assumed frequency. 
For symmetric bending, po is determined in such a way 


as to give zero slope at the base, and the quantity 
F, + Mywhy, (27) 


If the value of (27) is zero, the assumed 
Since 


is evaluated. 
frequency was the correct symmetric frequency. 
py) is known, the deflections at each station may be 
evaluated to give the corresponding mode shape. 
Similarly, for antisymmetric bending, fo is so deter- 
mined that the deflection is zero at the base, and the 
quantity 

M, - T,w" pp (28 
is evaluated. If (28) is zero, the assumed frequency 
was the correct antisymmetric frequency, and the cor- 
responding mode shape may then be evaluated. If 
neither (27) nor (28) is zero, the process is repeated with 
a new assumed value of the frequency. The correct 
frequency corresponding to each case may then be 
found graphically by plotting the values obtained for 
(27) and (28) as functions of the frequency. An alter- 
native method is to determine the quantities 7, or /,, so 
is to make (27) or (28) zero. This method is useful 
when the mass at the base is subject to variation due to 
different loading conditions of the aircraft, etc. 

The mass of a concentrated item at a point other than 
at the center may be introduced as a variable in the 
same manner as in the torsion case if rotary inertia is 
neglected. In this instance, all entries in columns 
4), (6), (10), and (15) at points inboard of the item 
will contain an additional term due to the concentrated 
mass. The mass and the tip angle then become two 
unknowns that may be simultaneously determined so 
as to satisfy the boundary conditions at the center. 

The torsional and bending computations are both 
greatly simplified if the stations are so selected that ( Ax) 
isconstant. For then, in the torsion case, column (9) 
becomes identical to (6), and column (10) is identically 
unity. In the bending case, column (6) may be re- 
placed by p,(Av/2); (7), by 2EI, +1; (8), by 
M,/EI,)(Ax/2)?; (9), by (EI, 41 + Fl,); (10), 
by F,( Ax/2)*; and (15), by A/,( Ax/2)?. Column (11) 
becomes identically 2 at all stations and may thus be 
dispensed with, while column (14) will no longer be 
needed. 

The insertion of two additional columns in the torsion 
calculations and one additional column in the bending 
provides a positive check at each station on the correct- 
ness of the results. To check the torsion case, a column 
is added in which the torque is evaluated by means of 
the expression 


T, = (GJ,)0,’ 
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TABLE 3 


Tabular Arrangement of Torsional Stiffness and Mass Moment 
of Inertia Data 










































































svat (LB.FT) x10” LB. SEC* FT: 
2 Gi, CIS ne, , Lr */ AXn 
o| 276 | .#6 | 204 | 336 |2.66667/-44= 
406 | 545 | 336 | 456 | 
2] 545 | 736 | 456 | 756 
S$ . 736 /,026 7.56 //.04 
— 
4] 74026 | 4472 | 672 | 20.60 
5S} 4472 2.268 /0.80 /2.44 
| A/RPLANE 
6 | 2.268 (244 -_— \Cenree NE 
CONCENTRATED NERT/A ” ; r < 
AT CENTRALING: 4200 LEFT, SEC. we 
—————-——_-——- - 4 
201 | eee ¢ oe 
/ 
> ———-—_— —-- . —-— -_ * 4 
| Gir, 
(44 séc*) | A y 
9 ——ES - 4 + 
GS*/C o Ps ‘ 
£8. Fr? 
.' i. —— | 
STATION No#O vf 2 3 re 5 6 
DISTANCE ) 
FROM Tips 0 32 6& 2% 28 92 
(/NCHES) | 


Figure 1 
Spanwise distribution of mass moment of inertia and 
torsional stiffness. 


Fic. 1. 


and a second where the torque is evaluated from 


_1+T7,0,) (29a) 


ro = —1 + w*( Ax,,/2)(I» _ 19, 
For the bending calculations, one column is added in 
which the shear is obtained from the relation* 


Pe = F,, —1 + w*( Ar, 2) (in = hy, —- 1 + Mnltn) (29b) 


In this way, any errors may be detected and corrected 
before proceeding further with the computation. 


(IV) NUMERICAL EXAMPLES 


(1) Torsional Frequencies and Modes with Variable 

Tip 

The semiwing is divided into six sections of equal 
width, Ax = 32 in., the points of division being num- 
bered from 0 to 6, with the point x = 0 at the wing tip 


* If (Ax) is constant, this may be written in the form 


Ax \3 Ax \3 Ax\4 
r{ 5) = F,, ‘ ( -) + w*( ) (un. thn _ i + Unltn) 
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Sample Calculation for Torsional Modes with Variable Tip 


TABLE 4 
























































Inertia; w?(Ax/2)? = 10 X 108 
a., * Gorn + Pia em 
(4) (5) (6) (7) (8) 
2) om. eat | G:() Jorn em [rast ane | 
o |~2°2-| 40396 <aaheg| cat |.20936 
! Sener | © 4164 200 951 | 50956 | 
2 |=S8278 -| 73/04. |2ES- 1.20/ | .74356 
3 | ser \.01844 |Z 1762 |/.03704 
4 | =$74¥ -|1.26520 Spamay| 2498 |/48260 
kei 99 2.25720 03624 | 3.740 \228 744} 
eg (a>. =03204 | 
9/37 Le 4.00007 I 
TABLE 5 


Calculation of Tip Inertia Necessary to Satisfy Centerline Condi- 


tions for Various Frequencies; Antisymmetric and Symmetric 



















































































Torsion 
AQT IGN a &% (ap) (ae 4 v0 Le, | Led, 
[ors] 25 | vos | steror) seer exter] 2725] 295 
}s53.| 50 | 2/0 segesl, ni 4 7, 704 ‘ase an 
[7s [700 | 020 |.aee, leans, [soteoen| 26 [ere 
Vos| 20.0 | ano | azees, | oer Se alan 
se| ea | soe | Bee, | bee] «S| 7 | 7 
|2/2| 00 | 93.60| 325/07, | nooeats,|n1a509 7, [Ow [O-6E 
|22s] 900 | 3700 | rroor, |no07dors.\-01s#3 7, pOoS |949-4 
2374 /00.0 | 42.00 200m, 007 Le Cees r, S662 | 5208 
|2594 /200 | 50.40 101934 Le * o06Nt T, “956/77, | | 20.95 | 2/.49 
bss 2000 | once | 75632, | _ase, [artes | ase | 230 
baal zea0| 22a |. Pe, | 2 | Zee | rac | rae 
Seza| 2400 |10080| , AES. | - 12001, |-5 72926. Z| -.086 | .06/ 
The tabulated values of J and GJ, listed in Table 3, 


are taken from Fig. 1. Since J is discontinuous at 
station 4, /,(m = 4) is taken as 6.72, while /, 4 1(” = 3) 
must be taken as 11.04 in Eq. (4a). Column (7) of 
Table 3 may now be written down and columns (5) and 
(8) computed for an assumed value of w 
w*(Ax/2)? = 10. The initial entry in column (4) is 1, 
and that in column (6) is found from the relation 


Ax T wo?( Ax/2)? 
(a)( )=-2 a ie 

2 GJo( Ax/2) 
resulting from the boundary condition at the tip. 
complete table may now be constructed by filling in 


= 75 sec.—! or 


— 0.002707, 


The 
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columns (4) and (6) with the aid of Eqs. (23) (cf. Table 
4). The boundary conditions at the base may now be 
applied. For antisymmetric torsion, the condition 


a = 0 


gives 


(Tq = 59.62 Ib. ft. sec.” 


For symmetric torsion, the condition 
GJy(0’) = Tw, 
results in 


(I,), = 61.4 Ib. ft. sec.’ 


By repeating the above process for other values of 

Table 5 may be constructed. 

Fig. 2 shows (/,), plotted as a function of w. The 
asymptote, corresponding to zero amplitude at the tip, 
may be located by interpolating for the point where 
the coefficient of J, vanishes. This point also marks 
the transition from the first mode to the second. From 
Table 5, it can be seen that, except for large values, 
the difference between (/,), and (J;), is negli- 
gible. 

The antisymmetric mode shapes corresponding to 
various values of (J,)4, obtained by resubstituting into 
the expressions for the deflections the appropriate 
values of this parameter, are shown in Fig. 3. To ob- 
tain the mode shape for the case when/,— @, it must 
be recalled that the tip deflection, while taken as unity 
for convenience, is actually arbitrary. Thus, all 
subsequent deflections contain this arbitrary quantity 
% as a factor. By allowing /, to become infinite and 
6) to approach zero in such a way that the torque 
wl My remains finite, it is seen that only the coefficients 


of [,;remain. These coefficients give the required mode 





shape. The symmetric mode shapes differ only 
slightly from the antisymmetric and are therefore not 
shown. 
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FIGURE S52 


Variation of the first antisymmetric torsional mode 
(J; is in Ib. ft. sec.?.) 


Fic, 3a. on of t! , 
shape with tip inertia. 
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FIGURE 3h 


Variation of the second antisymmetric torsional mode 
shape with tip inertia. 
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(2) Bending Frequencies and Modes with Fixed Tip Load 
and Variable Load at the Center 

The mass and bending stiffness distribution for this 
example is shown in Fig. 4. Ten equal divisions of 
width Ax = 1.9667 ft. are used, the basic data being ar- 
ranged in Table 6 in a manner similar to the data of the 
previous example. 

Columns (7), (9), and (13) of Table 7 are independent 
of the value of the frequency and may, therefore, be 
filled in at once, while columns (5) and (12) may be 
filled in for an assumed value. In the example, the 
quantity (w?/g)(Ax/2)* is used in place of the frequency 
itself, since Ax is constant and y is given in weight 
rather than in mass units. Table 7 is constructed with 
(w/g)(Ax/2)* = 0.3 X 10%. The initial entry in 
column (4) is unity and that in column (6) is a, where 
a = po( Ax/2) and fp is the arbitrary slope at the tip. 
Rotary inertia of the tip mass is neglected, and so the 
initial entries in columns (8) and (15) are zero. The 
initial value in column (10) is obtained from Eq. 
(25b), which is rewritten as 


Ax\3 Ax\*/ 2 
r(S) = ras ( >) (z) 


whence, 


Fo( Ax/2)? = 0.36183 X 108 


upon making use of the fact that w?/g (Ax/2)* = 0.3 X 
10°, mg = 1.186 X 10° and (Ax/2) = 0.98333. The 
table is then completed, using Eqs. (26) to compute the 
successive values in columns (4), (6), (8), (10), and (15). 
Each entry in these columns consists of two parts: 
one depending on the initial tip amplitude and the 
other on the initial tip angle. 

The values at the base in columns (4), (6), (10), and 
(15) are now used to determine the correctness of the 
assumed frequency. In the symmetric case, the condi- 
tion that the slope at the center be zero gives 





TABLE 6 


Tabular Arrangement of Bending Stiffness and Weight Distribu- 
tion 
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Figure 2 


Fic. 4. Spanwise distribution of mass and bending stiffness. 
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Sample Calculation for Bending Modes; 


(w?/g)(Ax/2)* = 


0.3 X 10%, 


TABLE 7 


Tip Weight.) 


(*Base Values Obtained by Including Rotary Inertia of 
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The second condition at this point—namely, 


F, = —mwh, 
gives 
mg = —12.3 lbs. 


In the antisymmetric case, the value of a which gives 


zero displacement at the center is 
a = —2.47675 


whence, from the second centerline condition, 


M, = Twp» 


it is found that 


I,g = 80.17 X 10° Ib. ft.? 


Repeating the process for other assumed frequencies, 
the variation of the symmetric frequency with mass 
at the base and of the antisymmetric frequency with 
the inertia at the base may be found. A convenient 
method of presenting these results graphically is as 
follows: 


tity 


F,(Ax/2)3/h, 


is evaluated for symmetric modes, and the quantity 


M,( Ax/2)?/p,( Ax/2) 


for antisymmetric (Tables 8 and 9). These ratios are 
plotted as functions of (w*/g)(Ax/2)* on an x-y coordt- 
natesystem. Inthe symmetric case, the above graphis 
intersected by a family of lines through the origin 
whose equation is given by 
y = [inng/(Ax/2) |x 

The points of intersection give the desired frequencies 
corresponding to varying values of mg. Similarly, the 
antisymmetric curves are cut by a series of lines defined 
by the equation 


y = [Lyg/( dx/2)3x 


In Figs. 5 and 6, the lines corresponding to the values | 


mg = 3 X 10,5 X 10° areshown. The radius of gyra- 
tion of the fuselage is taken to be 1.5 ft. so that J, = 


2.29 i». 
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TABLE 8 


Calculations for Construction of Bending Frequency Diagram; 


Symmetric Modes 
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Calculations for Construction of Bending Frequency Diagram; 
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From the figures and the data of Tables 8 and 9, it can 
be seen that, for7,g = 5,000 lbs., the first and second 
symmetric frequencies correspond to values of (w*/g) 
(Ax/2)4 & 10-3 0.006 and 0.21, while the first two 
antisymmetric frequencies correspond to (w*/g)(Ax/2)4 
X 10-* = 0.037 and 0.7. The symmetric frequencies 
and the first antisymmetric frequency are only slightly 
affected when M,g is reduced to 3,000 Ibs., while the 
value determining the second antisymmetric frequency 
is increased to 0.79. The actual frequencies are sum- 
marized as follows: 


For tig = 5,000 Ibs., 


84.9 sec.~! 
—] 


Ws, = 14.2 sec.—!, ws, = 


Wa, = 35.6 sec.—', wa, = 


For m,g = 3,000 lbs., 


150 sec. 


Ws, = 14.9 sec.— 3 sec.—! 


We, = 36.1 sec.—', ow, = 


Ws, = NT 
165 sec.—! 


frequencies for two values of i. (5a) 
and second modes. 


Fic. 6. 
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The above information may also be employed to show 
the effect of variable flexibility of the fuselage, such as 
would be encountered if the measured frequencies were 
obtained with the aircraft on partially inflated tires 
or air bags. Consider the symmetric case, and let 
w, be defined as the natural vertical frequency of the 
semifuselage itself on the flexible support. This quan- 
tity may be determined experimentaily by measuring 
the frequency of the entire aircraft on the support, thus 


! 
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LiIGURE © 
Location of the first and second antisymmetric bending 
frequencies for two values of J,. 
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Fic. 7. Change in the shape of the second antisymmetric 


bending mode due to rotary inertia of the tip mass 
determining the stiffness A,. The required frequency 
may then be found from the relation 


WwW” — K, Mp, 


The shear at the centerline must in this case be equal 
to the quantity 


Why + Kehr 


or, making use of Eq. (29), 


F "» = 


why [1 = (wy ‘w)?] 


Thus, the actual mass at the base may be replaced by 
an “‘effective’’ mass equal to 


mMy[1 — (w,/w)?] 


Since each assumed frequency w for the system de- 
termines a corresponding value for the mass at the base 
the result is equivalent to a fixed mass at the base with 
variable flexibility. 

The antisymmetric case may evidently be handled 
in an analogous manner. 


(V) Rotary INERTIA EFFECTS 


In general, the rotary inertia of elements of the wing 
itself is negligible. However particularly in higher 
modes of vibration, the inertia of concentrated items 
may have appreciable effect both on the mode shape 
and the resulting frequency. To illustrate this effect, 
the previous example was recomputed, taking into ac- 
count the rotary inertia of the tip weight. Only the 
second entry in each space is affected by the change. 
The initial entry in column (15) (Table 7) is deter- 
mined from the equation 
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M, = —I wp, 


(8) =e) QE) 


The radius of gyration of the tip weight about the 
bending axis is 1.225 ft., whence, 


Ax\? ___ Ley As? b 
M, —— = —1.779 = a X 10° 


The final values are shown in the last line of Table 7. 
Carrying out the correction for a second frequency 


or, 


corresponding to 


(aw? g)| Ax 2)4 —— 0.2 x 108 


determines the second curve indicated in Fig. 6. It can 


be scen that, for/,g = 11.25, the corresponding value of 


[(w?/g)( Ax /2)*] K 1073 

is decreased from 0.7 to 0.3, a reduction of 35 per cent 
in the frequency itself. It may also be noted that the 
inertia of the fuselage is less effective when the inertia 
of the tip mass is considered. In the symmetric case, 
there is also considerable reduction in the frequency of 
the second mode, although it is not so great as in the 
antisymmetric. 

The change in the shape of the second antisymmetric 
mode is shown in Fig. 7. Unlike the case when the 
rotary inertia was neglected, it is possible to obtain a 
mode shape with zero amplitude at the tip without an 
infinitely large tip mass. For this reason the deflec- 
tions were ‘‘normalized’’ to unity at the second inboard 


station rather than at the tip station. 
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A Note on the Velocity of Sound 


Morris Morduchow 
Polytechnic Institute of Brooklyn, Brooklyn, N.Y. 
June 25, 1949 


ci HE PURPOSE of this note is to show that the well-known dif- 
ferential expression a? = dp/dp for the velocity of sound can- 
not, as is sometimes attempted, be rigorously derived from the 
steady one-dimensional differential equations of nonviscous flow. 
It is also the purpose to show that these equations can, neverthe- 
less, be used to derive the expression a? = yp/p without the re- 
jection of any small-order terms. 

The basic differential equations for one-dimensional flow re- 
ferred to above and representing, respectively, conservation of 
mass, momentum, and energy are: 


pdu+udp =0 (1) 
pudu+dp=0 (2) 
udu+c,dT =0 (3) 


where u, p, p, and T denote velocity, pressure, density, and tem- 
perature, while c, is the specific heat at constant pressure. The 
general solution of Eqs. (1), (2), and (3) and the ideal gas law 
pb = pRT is the “‘trivial’’ solution u = constant, p = constant, 
p = constant, T = constant. One may, however, assume a dis- 
continuity present, separating two regions in each of which x, p, p, 
and T are constant, but with different values for the two regions. 
The derivative dp/dp is then zero everywhere except at the dis- 
continuity where its value is indeterminate. 

If the increments in u, p, and p are considered as finite though 
small but not as differentials, Eqs. (1), (2), and (8) are still valid, 
and it is possible by the above means to arrive at the result u? = 
Ap/Ap. It would, however, be incorrect now to let Ap — 0, to 
write a? = lim (Ap/Ap) = dp/dp, and to interpret dp/dp as a 
derivative, since in the limit of the equations of motion reduce to 
the differential forms (1), (2), and (3), and the false implication 
would have been made that the derivative dp/dp is determinate at 
the discontinuity. 

It may be interesting to observe that attempts to derive 
formulas by ignoring the indeterminateness of the derivative at 
the discontinuity will lead automatically to mathematical incon- 
sistencies. For example, suppose that from Eqs. (1) and (2), one 
Writes 


u? = dp/dp (4) 
while from Eqs. (2) and (3) and the gas law one obtains 
pb = kp” (5) 
where & is an arbitrary constant and y is the ratio c,/c, of specific 
heats. Then considering (incorrectly) Eq. (4) as a derivative, 
it would follow from Eqs. (4) and (5) that 
yb/p (6) 
If, however, to be consistent, expression (6) is now considered as 
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4? = 





differentiable for substitution into Eq. (4), giving 
(1/2) yd(k/p) 
substitution into Eq. (2) would lead to 


udu = 


p = cpv/(y+2) (7) 
where c is an arbitrary constant. The change of state [Eq. (7)] 
evidently contradicts Eq. (5). 

In spite of such inherent inconsistencies contained in treating 
Eqs. (1), (2), and (3) as if derivatives such as dp/dp existed, the 
above type of derivation of expressions (4)—(6) is sometimes seen 
in the literature. An tlegant way of rigorously deriving Eq. (6) 
from the (physically correct) differential equations [Eqs. (1), (2), 
and (3)] is to integrate these equations. Thus, these equations 


become: 
pill; = potle (la) 
pitts? + pi = potte® + pr (2a) 
(%:2/2) + ¢pT = (t2?/2) + ¢yT2 (3a) 


where the subscripts 1 and 2 refer to values at two different points 
in the field. From Eq. (la), (2a), and (8a) and the gas law, it is 
found that the pressure ratio p/p, must have either a value of 
identically unity (‘‘trivial’”’ solution) or a value (Rankine-Hugon- 


iot) of 
p_ 2 (=) _y-1 o- 
pa oyt+l\ pb y+1 


Eq. (8) represents a wave of finite discontinuity (‘shock 
wave”). Defining a sound wave now as an infinitesimally weak 
shock wave—i.e., as one for which ~./p; — 1—it follows from 
Eq. (8) that for a sound wave 





pilly?/py > y or m2 > ypi/p1 


Thus expression (6) for the velocity of sound is derived in a neat 
rigorous manner from the straightforward differential equations 
of motion, without the rejection of any terms in these equations. 
The differential expression (4), however, cannot be so derived. 
A proper derivation of Eq. (4), based on the equations of unsteady 
flow, is given, for example, in reference 1. 
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Lift and Drag of Rectangular Wings of Finite 
Aspect Ratio in Accelerated, Transonic 
Flight 


C. S. Gardner and H. F. Ludloff 
New York University, University Heights, New York, N.Y. 
July 8, 1949 


t N A PAPER PRESENTED at the I.A.S. Seventeenth Annual Meet- 
ing in January, 1949 (Preprint No. 186), the authors investi- 
gated the influence of acceleration on lift and drag of a two- 
dimensional airfoil moving at supersonic or transonic speed. In 
particular, it was shown how the aerodynamic characteristics re- 
main finite and vary continuously when the flight speed U + bt 
passes through the sound speed a. The purpose of the present 
note is to indicate how the method and results can be extended to 
wings of finite aspect ratio. In the latter case, the applicability 
of the linearized theory in the transonic range is more legitimate 
than in the case of infinite aspect ratios. 
The drag can be derived by use of retarded potentials, which 
can be brought into the form of Possio salutions: 


-f fy a dtd 
' Vit —*- (1/e%%-)* 
R 
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where o@ is given by 


o(f,7) = —(1/r)¢y —(1 w)f’ [é + Ur + 


f(x) being the airfoil profile form. 

One can show that, for rectangular wings of finite aspect ratio, 
the domain of integration in the x t-plane, R, is given by the 
shaded area in Fig. 1, which is bounded by the two sound speed 
lines, passing through any point 7 on the airfoil surface, by the 
parabolic path of the leading edge and by a hyperbola whose 


h = span length of wing) from P. 


vertex is at a distance }b/a (b = 
(See Fig. 1. 


produced by the airfoil is small compared to the dynamic pressure 

at the leading edge—in other words, as long as C, < 1 or, if we 

average the pressure distribution over the chord length, sa <q 4% 

Taking a thin supersonic airfoil of 5 per cent thickness ratio, we 
have, at JJ = 1, for B = 0.02 

- ~ 0.25 for A.R 

0.2 for A.R.7 

0.14 for A.R. 2 


Q 


Therefore, the above results might be considered as reasonably 
valid, even at J = 1; the smaller the A.R. is, the better becomes 


the approximation. 
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Effect of Disposable Load on Range of 
Turbojet Airplanes 


W. F. Savage 
Assistant Professor, Department of Mechanical Engineering, 
University of Kentucky, Lexington, Ky. 


August 1, 1949 


N CONSTANT ALTITUDE CRUISING FLIGHT of turbojet-propelled 
I military aircraft, analytical equations can be derived for the 
range and average operating speed when a bomb or cargo load is 
dropped at the midpoint of the flight. 

As outlined by Ashkenas,! a formula can be derived for the 
maximum range of turbojet aircraft. By considering flight at 
constant C, Cp, which would be possible for most aircraft at 
altitudes equal to, and below, combat ceiling, the formula for the 
maximum range of an airplane for a ferry mission would be: 

Rimiles) = 39.6 - = : (w,'/? — W,'/2) 

Cp (eS) “* ¢ 

wing area, o is the altitude density ratio, c is the 
fuel consumption, W, is the take-off gross weight, 
and IW’, is the landing gross weight. The thrust specific fuel con- 
sumption is assumed essentially constant, and the entire weight 
difference between the take-off and landing gross weights is con- 


where S is the 


thrust specifi 


sidered due to fuel only. 

The above equation is used to develop the range expression for 
the condition where a bomb or cargo load is dropped at the range 
midpoint. If WW’, is the gross weight of the airplane before drop- 
ping bombs or cargo and if WW’, is the gross weight of the airplane 


after the drop, the maximum range equation can be written as: 


~ 1 
; C, ”’ l ee ‘ 
R(miles) = 2(39.6) — = - (W,'/? = W,'*) 
c D ig>)‘* € 
or 
oo 
' Ge"? l l . : 
R(miles) = 2(39.6) — —~ - (W,'/? — Ww,’ *) 
Cp (a5) ¢ 
since the outgoing range leg is equal to the return range leg. The 
substitution that W, = W, + Wp can also be made, where Wp 


is the cargo or bomb weight and expressions for W, and W, can 
be written in terms of W;, Wy, and Wp. 


W, = [(N2 + Wp)/2N]?2, W, = [(N? — Wp)/2N]? 
where VN = W,'/? + W,'7?. 
Since the take-off gross weight is the sum of the disposable 
weight, the landing weight and the fuel weight, or (W; = Wp + 
Wy + Wer), the disposable weight can be eliminated from the 


range equation by inclusion of the fuel weight, Wr. Thus, 


w. = | 22+ (Ws — Ws — We) P 
2 2N : 
~ tea W, — Wr)? 
a. aN 


If the equation for either W, or Wy, is substituted in the range 
equation with bombs or cargo dropped, the resulting equation 
can be simplified by introducing the fuel weight-take-off weight 
ratio, Ay; the landing weight-take-off weight ratio, K; and the 
take-off wing loading, m. The final range equation for the condi- 
tion where bombs or cargo are dropped at the range midpoint will 


Cr (") | Kr 
Cp o c 1 a K' . 


This equation does not consider any expenditure of fuel for warm 
up, take-off, or climb, but this can be easily allowed for by con 
sidering Il’; as the gross weight when cruise begins at altitude. It 
will be noted in the above equation that the disposable load does 


be obtained 


Ri miles) = 39.6 


Not appear, but the equation is only applicable to the condition 
where the disposable load is dropped at the range midpoint and 
the fuel load is based on such an operatidn. 


READERS’ 





FORUM 37 





The average operating speed for the condition where a dispos- 
able load is dropped at the range midpoint can be written as: 
[Kre/(l + K”)] 


[l(2Kat K+ Ky 
nj — pares ae te 
K\ 2K"? +2— Kr 


where JV; is the initial cruising speed at take-off gross weight. 

For the case of constant speed cruising flight, the range formula 
for the condition where a disposable load is dropped at the range 
midpoint can be derived. 





Vave- = 4V; 





VC W, 
R(miles) = 2— 2 In — 
€ Up Vu z 
or 
VC, VW 
R(miles) = 2——* in — 
c Cp VW S 


Since Wy = W, — Wp, it is possible to solve for Wz in terms of 
Wi, Wy, and Wp. Thus, 


W, = (Wp + V Wp? + 4W/W))/2 


and 
VC 2W. 
Riniles) = 3 —-— 2 — 
cCp  Wo+V Wp? +4WyW, 


In the constant speed condition, the lift and drag coefficients 
will vary as the gross weight changes but it has been found by the 
author that the use of an average value of these coefficients based 
on landing and take-off gross weights is sufficiently accurate for 
most range calculations for constant speed cruise. 
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An Empirical Study of Low Aspect Ratio Lifting 
Surface with Particular Regard to Planing 
Craft 


N. Scholz 

Dr.-Ing., Institute of Fluid Mechanics, Technical University, 
Braunschweig, Germany 

August 1, 1949 


I AN ARTICLE by F. W. S. Locke, Jr. (Journal of the Aero- 
nautical Sciences, Vol. 16, No. 3, p. 184, March, 1949), the 
author attempts to solve the problem of the lifting surface of 
small aspect ratio purely empirically, using given measurements. 
But when giving empirical formulas, one should not overlook the 
theoretical results obtained with great difficulty. One should 
respect them in the structure of the empirical equations. 

The theory of the lifting surface has been developed as a linear 
theory, and, within the scope of its validity, it agrees well with 
Compared with a theory of higher order, 
it gives at least the lift and pitching moment slope for zero inci- 
dence. It does not contain the influence of the finite wing thick- 
ness on the total lift. In setting up an empirical equation for the 
lift slope the following points should, therefore, be taken into 


experimental results. 


consideration: 

(1) For a — 0, the formula should reduce to the linear lifting 
surface theory. 

2) It should be set up for the flat plate (infinitely thin). 
The finite thickness and the shape of the profile should then be 
accounted for by means of an additional factor. 

The lift coefficient should therefore be developed in powers of 
sin a, but one needs only consider the first two terms. For an 
extremely thin flat plate, one obtains, then: 
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Cy, = K, sin a + Ke sin’? a (1) 
K, and K, are dependent on the aspect ratio, as well as on the 
plan form of the lifting surface. Their limiting values are: 
A—0: K,=0; K; = 2 
Ao: K, = 2x; Ke=0 
K, can be obtained from the linear theory of the lifting surface 
for every aspect ratio and plan form. For an elliptical plan form, 
the following formula can be given as I have shown in a recent 
paper:! 
K, = 29A/(2 + Wa? + 4) (2) 
For A — 0, there is: 
K, = (#/2)A (3) 


Eq. (3) is valid for all plan forms, as has also been given by 
Jones.2, With good approximation it can be used for aspect ratios 
A <1. Bollay’ calculated the value of Kz = 2forA —0. This 
result is in good agreement with measurements. For aspect 
ratios A > 3, the curves Cz(@) are always linear, and therefore 
K, = O for A > 3. 

It is well known from measurements at large aspect ratios that 
the influence of the finite thickness of the wing and the profile 
shape on the lift is extremely small as long as the thickness ratio 
remains within normal limits. For aspect ratios A < 3, such an 
influence is, however, present, as has been shown by different 
measurements. Compared with flat plates, the lift slope of air- 
foils for a = 0 is smaller. The influence of the finite thickness 
would have to be accounted for by means of a correction factor 
for the constants K, and Ko. 

However, there are still other influences on the lift. The flow 
around the wing tips is of great importance. This flow, on the 
other hand, is mainly dependent on the shape of the side edges 
and on the wing plan form and profile. As the situation is not 
at all clear, only extensive systematic pressure distribution 
measurements may solve this problem. That there is such an 
influence is, however, clearly shown by the measurements pub- 
lished by Winter‘ and Zimmermann.’ However, not enough 
reliable measurements have been made to consider these effects in 


our lift formula. 
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Base Pressures at Supersonic Velocities 


Ingénieur Général Gabeaud 
Centre d’Essais de Mécanique des Fluides, Paris, France 
July 9, 1949 


7 THE Marcu, 1949, IssuE, the JOURNAL OF THE AERONAUTICAL 
SCIENCES published, under the title ‘‘Base Pressures at 
Supersonic Velocities,’ a well-documented paper of Freeman K. 
Hill and Ralph A. Alpher. I read on page 159, that “there has 
not been published any satisfactory theory to predict the values 
of the base pressure. However, in 1932, von K&rman and 
Moore developed a semiempirical formula that gave the base 
pressure coefficients as 
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VOW KARMAN & MOORE FORMULA 
AVERAGE APL-JHU FREE FLignT 


BASE PRESSURE COEFFICIENT, 


MACH HUMBIR. 


Fic. 1. 


b/po = {1 — [(y — 1)/2]AM24V OY 


where p is the absolute value of the pressure at the base, pp is free 
stream static pressure, y is the ratio of specific heats, and \ is an 
arbitrary constant to which on the basis of low subsonic measure- 
ments von K4rm4n and Moore assigned the value 0.2.” 

I wish to point out, for what use it may have, that I published 
in 1931 a turbulence theory of residual resistance.! This theory 
is, moreover, known to Dr. von Karm4n, who referred to it in the 
paper he presented October 2, 1935, at the Volta Congress;? at 
the same conference (p. 281) Margoulis pointed out that I had 
subsequently made an “improvement”’ in my theory. Returning 
to this subject, in 1936 I gave new developments and presented 
a curve giving unit base pressure as a function of speed. If the 
law of variation of pressure coefficient as a function of Mach 
Number is inferred from it, there is obtained the curve repre- 
sented by a heavy line on Fig. 1, where I have reproduced the two 
curves (free flight and von Karman and Moore formula) taken 
from Fig. 10 of Hill and Alpher’s paper. It may be maintained 
that my theoretical curve of 1936, when superimposed on the 
average experimental curve of free flight, presents a significant 
correlation, notably in the region (0.9 < M < 1.7) where disagree- 
ment is particularly obvious between the free-flight results and 
those of the semiempirical formula of von Karman and Moore. 

I wish to add that, since 1936, I have again added refinements 
to my theoretical outline; the results, not published, coincide still 
further with those of free flight. 


REFERENCES 
1 Comptes-rendus del’ Academic des Sciences, Vol. 192, p. 1630, 1931 
2 Reale Accademia d'Italia, Fondazione Alessandro Volta, Convegn 
Scienze Fisiche, Matematiche e Naturali, 1935. Le Alte Velocita in Aviasione 
pp. 270, 271, 281; Rome, 1936 
3 ‘‘Recherches sur la résistance de 1’Air,”’ (Memorial de l'Artillerie Fran 
¢aise, Tome XV, 4°™* fascicule de 1936, pp. 1290-1300) 


-_-— _ = _ — 
Propeller Flutter 
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July 21, 1949 


foes AUTHORS would like to clarify the simplifications made 1 
the body of their paper! regarding the mode shapes to be 
employed in a propeller flutter analysis. It is to be recalled that 


* Presently employed as Vibration Engineer at the Boeing Airplan¢ 
Company. 
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READERS’ 


the normal modes there derived were based on a clamped root 
section in the hub. The emphasis to be placed on appropriate 
displacement forms in a propeller flutter study stems from a wide 
background of experience in attempting to represent the vibra- 
tory characteristics of a complex beam by as few normal modes as 
possible. Small errors in this representation have large effects 
on the important energy forms and often produce results that are 
subject to question. 

In generalizing tHe analysis presented in the paper, it is desir- 
able to consider the entire vibrating system of the propeller and 
its power plant. Because the blade is twisted, the normally un- 
coupled axial and torsional motions of the propeller shaft result 
in coupled modes in the blade itself. The resulting “system 
modes’’ can be classified in the following manner. 

Class I—Torquewise angular displacements of the propeller 
hub are accompanied by linear axial motions due to blade twist. 
A simplification of this class forms the basis on which the original 
derivations were made. 

Class [IJ—Whirling modes that result when work is done on the 
propeller by the rotating forces within the engine. Because of 
blade twist, angular motions of the shaft axis are coupled with 
lateral motions of this axis so that the centerline describes a 
hyperboloid of revolution. 

Class IJJ—Reactionless modes that result from certain phase 
relations between blades at multiples of propeller rotational 
speeds. At the center of rotaticn these result in complete can- 
cellation of vibratory reactions. 

Extensive subdivision can be made in each of these classes for 
particular conditions. An exact treatment of any type of system 
modes would be accomplished by combining the engine and pro- 
peller components of the system through an impedance concept 
in order to make engine and propeller motions and reactions con- 
formable at the blade root. Since this treatment would require 
the consideration of mathematical complexities not directly 
associated with the essential features of propeller flutter, it was 
not included in the paper. 


REFERENCE 


‘Turner, M. J., and Duke, James B., Propeller Flutter, Journal of the 
Aeronautical Sciences, Vol. 16, No. 6, pp. 323-336, June, 1949. 
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Compressibility Correction for Axially 
Symmetric Bodies 


E. V. Laitone 


Associate Professor, University of California, Berkeley, Calif. 
Also Consulting Aerodynamics Engineer, Hughes Aircraft 
Company 

July 18,1949 


Q. (14) OF REFERENCE I presents an integral that is non- 
existent and therefore cannot be considered proof that there 
is no compressibility correction on an axially symmetric body. 
The correct expression for the surface pressure coefficient on a 
symmetrical body of revolution with a local cross-sectional area of 
S(x) and a total length L is given in reference 2 for linearized sub- 
sonic flow as 
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where 8 = V 1 — M2andr™ = dr/dx. 

This shows that the compressibility effects depend on the actual 
location on a given body, as well as the shape of the body as de- 
fined by the second derivative of the cross-sectional area. For 
example, at the midpoint of any symmetrical body of revolution 


1+ (108 $14 106. -- 
1). 


L? S® 
48 §(2) 
As seen in Fig. 2 of reference 2, the compressibility effects are 
negligible for a prolate spheroid less than 5 per cent thick pro- 
vided the Mach Number remains less than 0.8. 
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Erratum 


W. E. Strohmeyer and D. R. Gero 

North American Aviation, Inc., Los Angeles, Calif., and Goodyear 
Aircraft Corporation, Akron, Ohio, respectively 

September 13, 1949 


TT FOLLOWING ERROR appeared in our Readers’ Forum item 
“An Approach to Supersonic Airfoil Theory” (JouRNAL OF 
THE AKRONAUTICAL SCIENCES, Vol. 16, No. 8, pp. 506-507, Au- 
gust, 1949). 
The first equation in the right-hand column at the top of page 
507 should read: 


A P.E. = (7 + 4 


Jae Av = Udi 
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of the Libraries of the 


Institute of the Aeronautical Sciences 





The services of the Libraries are available ro all 


bers of the Institute, to Corporate Members, to advertisers 
Review and Agro- 


mem- 
in the AERONAUTICAL ENGINEERING 
NAUTICAL ENGINEERING Catatoa, and, under usual library 
Four specialized services are 


to the public. 


limitations, 
available 


The Paul Kollsman Lending Library 


This lending library service makes available, without 
charge, the latest and more important aeronautical books. 

Members may request the loan of any aeronautical or 
technical book they wish to borrow. Through an exchange 
agreement with the Engineering Societies Library, any 
book on general engineering may be borrowed from its great 
collection of over 160,000 volumes. 

A photostating service is available at usual library rates. 

Applications for membership in the library and further 


information will be sent on request 


The W. A. M. Burden Reference Library 


This reference library contains over 12,000 acronautical 
books, magazines, pamphlets, and reports gathered from 
world-wide sources and is one of the most complete aero- 
nautical libraries in the world. Material from this library 
is not available for loan but may be used for reference pur- 


poses. 


The Pacific Aeronautical Library 


7660 Beverly Boulevard 
Los Angeles 36, California 


Established in cooperation with the aircraft companies 
the library serves. The leading aircraft companies in or 
near Los Angeles participate in its support and operation. 

This service library for aeronautical research is available 
to the public for reading privileges. Source material in- 
cludes aerodynamic and structural research reports, as well 


as books on drafting, production methods, history, and al- 
licd sciences. It furnishes books, periodicals, and pamphlet 
material to the participating aircraft companies to supple- 
ment their engineering libraries. 


Technical Information Service 


This service has experienced personnel under the super- 
vision of trained aeronautical engineers to compile any in- 
formation desired. The services range from listing special- 
ized reference books to the preparation of exhaustive 
bibliographies, digesting of reports, and general surveys of 
any acronautical subject. Some of the available services 


are: 


Bibliographies on any aeronautical subject. 

Reports on any aeronautical subject. 

Digests of aeronautical books, papers, periodicals, and 
ences. 

Translations. 

Engineering investigations of special acronautical subjects. 

Biographies of individuals engaged in aeronautics 

Photostats of any aeronautical or general engineering 
material. 

Microfilms made on special order. 

Photographs made from the Institute's photographic collec 
tion. 

Drawings and tracings made. 


In addition to the services mentioned any commission 
which comes within the scope of the Service will be ac- 
cepted. Special arrangements may be made for work re 
quiring several weeks or months. 

Translators are available for accurate transcriptions of all 
foreign language data. Translations are carefully edited 
by trained engineers. 

Reproductions of any matcrial in the collections of 
the Institute may be ordered at standard photostat rates. 
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2 EAST 64th STREET, NEW YORK 21, N.Y. 
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